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LIQUID JET IMPINGEMENT NORMAL TO A DISK IN ZERO GRAVITY 


by Thomas L. La bus 
Lewis Research Center 

S UMMARY 

An experimental and analytical investigation was conducted to determine the free- 
surface shapes of circular liquid jets impinging normal to sharp-edge disks under zero- 
gravity conditions. An order- of- magnitude analysis was conducted that indicated re- 
gions where viscous forces were not significant in the computation of free- surface 
shapes. The demarcation between the viscous and inviscid regions was found to depend 
on the flow Reynolds number and the ratio of the jet and disk radii. 

Experiments conducted under zero-gravity conditions yielded three distinct flow 
patterns. These flow patterns were defined as surface tension flow, transition flow, 
and inertial flow. The flow regions were classified in terms of the relative effects of 
surface tension and inertial forces. The transition between regions was correlated with 
the system Weber number and the ratio of the jet and disk radii. 

A zero-gravity inviscid analysis was performed in which the governing equations 
and the boundary conditions in the physical plane were transformed into an inverse 
plane. In the inverse plane, the stream function and the velocity potential became the 
coordinates. This removed the prime difficulty in free- surface problems, that of 
having to guess at the true position of the free surface. The governing equations were 
nonlinear in the inverse plane, thus requiring a numerical solution in which sets of 
nonlinear algebraic equations were solved simultaneously. Comparisons between ex- 
periment and numerical computations were made for the infinite- and finite-plate cases 
with the result that good agreement for the free-surface shapes was obtained. 


INTRO DUCTION 

A knowledge of the dynamics of free liquid jets is required for the solution of a 
variety of problems associated with fluid flow within propellant tanks under low gravita- 

* Submitted as a thesis in partial fulfillment of the requirements for the degree 
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tional conditions. In particular, an understanding of the liquid-jet impact process, 
such as occurs when a liquid impinges upon baffles or tank walls diming an inflow or 
reorientation maneuver, is required to predict liquid- propellant location, heat transfer 
rates, and pressure distributions. Liquid jet impingement also has direct applicability 
to spacecraft fire safety for cases in which water jets are used to extinguish fires 
under low-gravity conditions. To predict delivery flow rates requires the accurate 
prediction of flow surface coverage as a function of jet momentum flux. 

There generally are three chief obstacles that have in the past prevented the solving 
of steady-state liquid jet - solid interaction problems. The major obstacle is the 
presence of the free surface. To apply numerical techniques to the solution of free-jet 
problems, it is necessary to define the area over which the computations are made by 
means of boundaries determined by the free liquid surfaces. Unfortunately, the location 
of the free surface is one of the solutions sought, so various complex techniques must 
be devised to circumvent this situation. Furthermore, analytical techniques are re- 
stricted solely to two-dimensional problems, whether a free surface exists or not. The 
second obstacle is gravity. Liquid jets in air (free jets), unlike liquid- into- liquid jets 
and gas-into-gas jets (submerged jets), are affected significantly by gravitational forces. 
The free- surface shape and velocity profiles depend on both the magnitude and orienta- 
tion of gravity. Adding gravity necessarily complicates a model either through the 
governing equations or through the boundary conditions. Neglecting gravity in the model 
makes questionable the comparison of the theory with normal- gravity experimental 
data. The final obstacle is surface tension, an effect that has generally been neglected 
in almost all studies on free jets. Adding surface tension into a model leads to non- 
linear free- surface boundary conditions. 

This report presents the results of an experimental and analytical study conducted 
at the NASA Lewis Research Center concerning zero-gravity isothermal liquid-jet im- 
pingement. An axisymmetric liquid jet was impinged normally onto a sharp-edge disk 
under conditions in which both inertial and surface tension forces are of importance. 

The experimental free-surface shapes were correlated with known system parameters. 
An analytical model was formulated, and the free-surface shapes and streamlines were 
calculated for a number of discrete cases. 


LITERATU RE SURVEY 
EXPERIMENTAL STUDIES 

Very few experimental studies have been conducted to examine free jets impinging 
on solid surfaces. No work has been done where the major concern was either the shape 
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of the free surface or the measurement of velocity profiles within the jet. Also, only 
one experiment has used a two-dimensional jet. A two-dimensional jet is one in which 
the flow emanates from a rectangular slot so that the width of the jet is very large rela- 
tive to its thickness. Schach (ref. 1) measured the pressure distributions and analyti- 
cally calculated the free-surface shapes and velocity distributions for jets impinging 
onto flat panels at various impingement inclinations relative to the direction of flow. 

The jet employed had dimensions of 21 millimeters by 115 millimeters. According to 
Schach, the jet diverged spatially after impinging on the panel, and thus it can only be 
considered as truly two dimensional close to the centerline. 

An excellent account of an elaborate experimental apparatus for obtaining a quies- 
cent circular water jet in normal gravity is given by Donnelly and Glaberson (ref. 2). 
Their major concern was jet stability under imposed audiofrequency disturbances and, 
therefore, the impingement phenomenon was not directly observed. Rupe (ref. 3) and 
Stephens (ref. 4) experimentally measured the pressure distribution caused by circular 
jets striking solid surfaces in normal gravity. However, neither Rupe nor Stephens 
measured the free-surface shape or discussed any instabilities that occurred. 

In nearly all flows where a circular liquid jet strikes a large flat surface in normal 
gravity, the impinging liquid jet moves radially outward from the stagnation point until 
a certain radial distance is reached, whereupon an instability known as a circular hy- 
draulic jump occurs. The jump is characterized by an abrupt increase in the liquid 
depth and turbulent fluid motion. Koloseus and Ahmad (ref. 5) were concerned solely 
with predicting the behavior of the circular hydraulic jump. A water jet impinging on a 
flat plate of epoxy material was used in these experiments. The circular hydraulic jump 
was the subject of a very complete study conducted by Nirapathdongporn (ref. 6), whose 
report contains an excellent description of various devices for measuring jet shapes and 
jet diameters. 

All of these studies deal with normal- gravity liquid jet - solid impingement. There 
have been no experimental studies on the impingement of liquid jets under zero-gravity 
conditions. 


ANALYTICAL STUDIES 
Steady Two-Dimensional Potential Flow 

A number of papers and books have presented analytical studies of steady-state, 
two-dimensional free jets impinging on a variety of surfaces. The majority of these 
studies were concerned with irrotational, incompressible, inviscid flow in which the 
effects of gravity and surface tension were neglected. One of the major attractions of 
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this approach to the problem is that it can be handled by using complex potential theory 
and, therefore, can be treated analytically. 

A two-dimensional jet striking an infinte flat surface at various angles was exam- 
ined by Batchelor (ref. 7), who solved for the limiting stream thickness as a function 
of flow impingement angle and jet diameter. However, no attempt was made to predict 
free-streamline shapes. Schach (ref. 1) treated the impingement as a function of angle 
by using Prandtl's hodograph method and obtained the equations for free-surface shapes, 
flow distribution, and pressure distribution for the case of impingement on an infinitely 
wide plate. Kochin, et al. (ref. 8) also examined the impingement of a two-dimensional 
jet obliquely to an infinite flat plate and discussed the case of impingement on a plate 
of finite length. The equation of the free surface for the case of a two-dimensional jet 
striking a flat surface at right angles is presented by Milne-Thomson (ref. 9), who also 
solved for the velocity components within the jet. An excellent discussion of the tech- 
niques for handling two-dimensional free- jet problems is presented by Gurevich 
(ref. 10). Some of the two-dimensional flows examined by Gurevich include flow around 
a finite wedge, flow perpehdicular to a finite plate, flow oblique to an infinite flat plate, 
and flows where a variety of solid objects are positioned adjacent to one wall or between 
two walls. Chang and Conly (ref. 11) analyzed the two-dimensional flow of a jet inter- 
acting with a number of flat segments at angles to one another. The results include 
flow turning angles but not free-surface shapes or velocity profiles. The irrotational 
flow pattern of a free jet discharging from a slot and flowing past a wedge was analyzed 
by Arbhabhirama (ref. 12). 

All of these papers and books are concerned with analytical techniques for obtaining 
solutions. The area of steady two-dimensional potential flow represents the most com- 
plete area of research in the field of jet impingement. 


Steady Axisymmetric Viscous Flow 

Watson (ref. 13) has analytically investigated free-jet impingement for the case of 
large Reynolds numbers where the viscous forces are confined to a thin boundary layer 
adjacent to the plate. A similarity solution was obtained for both the two-dimensional 
and axisymmetric velocity profiles and free-surface shapes for the case of normal im- 
pingement. As mentioned by Watson, the similarity solution can only be expected to be 
valid when the radial distance is sufficiently large for the incident jet to have lost its 
influence. The effects of gravity and surface tension were neglected in the analysis. 
Watson solved for the radial position of the circular hydraulic jump. 
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NUMERICAL STUDIES 


Steady Two-Dimensional Potential Flow 

When the shape of the solid upon which the jet impinges becomes complex, numeri- 
cal techniques for the solution of free- jet problems have to be applied. Jeppson 
(ref. 14) presents an excellent article in this regard. Jeppson used the stream function 
and the velocity potential as the independent variables and the coordinates as the de- 
pendent variables. A similar inversion approach has been used to solve a variety of 
fluid dynamics problems, as shown in references 15 to 19, and is mainly attributable to 
Thom and Apelt (ref. 15). Using this technique, Jeppson was able to circumvent the 
problem of working in the physical plane and having to guess at the true position of the 
free surface. The latter iterative approach was used in references 20 to 26 with limited 
success. Jeppson solved the problem of two-dimensional flow over a wedge and is the 
only one to have attempted numerical solutions of this problem. Lastly, Chan (ref. 27) 
applied the finite- element method to a number of free-surface flow problems, including 
flow from a circular orifice. 


Steady Axisymmetric Potential Flow 

Potential axisymmetric flow problems cannot be solved by using the powerful tool 
of complex analysis. Extending the mathematical treatment to axisymmetric and three- 
dimensional flows has so far proved intractable (ref. 9). For this reason, only numer- 
ical solutions can be attempted for problems of this nature. 

LeClerc (ref. 28) studied the impingement of an axially symmetric liquid jet normal 
to a flat surface. The shape of the free surface was found by using an electrical analogy. 
This method thus fixed the position of the free surface and enabled the author to apply 
standard finite- differencing methods and to employ Southwell's relaxation technique in 
order to solve Laplace’s equation. Jeppson applied his inversion technique to find the 
flow pattern and the free-surface shape for the case of axisymmetric flow past a variety 
of bodies of revolution, including cones. Jeppson applied his own technique to the solu- 
tion of a jet of inviscid, incompressible fluid issuing from a nozzle into the free atmos- 
phere. He shows how this method can be extended to a variety of other problems. 

Schach (ref. 29) used a semianalytical technique based on Trefftz's approximate method 
to find the shape of an axisymmetric free jet impinging normally on a plate. Also pre- 
sented in Schach's article was the pressure distribution on the plate, which was calcu- 
lated from the velocity distribution by using Bernoulli's equation. Young, et al. (ref. 30) 
and Brunauer (ref. 31) determined the flow pattern past two disks immersed in axisym- 
metric flow. 
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Both Young and Brunauer solved Laplace's equation in the physical plane. No 
analyses have been conducted for the case in which an inviscid free jet impinges on a 
plate of finite diameter. 

References 28 to 31 encompass all the known solutions with regard to axisymmetric 
jet impingement. References 18, 21, 22, 25, and 27 deal specifically with numerical 
methods applied to free- surface problems in which no impingement occurs. Jeppson 
(ref. 18) used an inverse formulation; the others worked in the physical plane. 


Unsteady Two-Dimensional and Axisymmetric Potential Flow 

Huang, Hammit, and Yang (refs. 32 and 33) have investigated unsteady flows and 
considered the impact phenomena for both two-dimensional and axisymmetric jets. The 
major interest in these articles was in obtaining the initial pressure distribution due to 
liquid impact. 


Steady Potential Flow Including Gravitational Effects 

Adding gravity in numerical analyses for potential flows causes no serious formula- 
tion problem for either the two-dimensional or the axisymmetric case, since its effect 
enters only through the free-surface boundary conditions and not the governing equa- 
tions. Jeppson (ref. 14) included gravity in his analysis of impingement on a two- 
dimensional wedge. Moayeri and Strelkoff (ref. 19), Southwell and Vais ey (ref. 22), 
and Chan (ref. 27) all considered the effect of gravity in dealing with steady, potential, 
free-surface problems in which no impingement occurs. 


Steady Potential Flow Including Surface Tension 

Zhukovskii (ref. 34) has indicated how to include the effects of surface tension. He 
examined a two-dimensional problem by using complex analysis, but his method is not 
extendable to either axisymmetric or three-dimensional flows. 


Steady Three-Dimensional Potential Flow 

Until very recently, little had been accomplished in the area of three-dimensional 
potential flow with a free surface, much less including impingement. Davis and Jeppson 
(ref. 35) developed a computer program to solve free-surface problems of this type by 
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the inverse method. Michelson (refs. 36 and 37) also examined jets under these condi- 
tions. He treated the case of an axisymmetric jet impinging obliquely on a flat surface 
and analytically showed the occurrence of wedge-shaped dry zones when the impingement 
angle was less than a critical value. Free-surface shapes are not obtainable by 
Michelson's method. 


ORDER - O F- MAGNITUD E ANALYSI S 
GENERAL FORMULATION 

The problem under consideration is the viscous flow of a circular liquid jet as it 
impinges normally to a flat plate, as shown in figure 1. The objective is to determine 
the free-surface shape of the impinging liquid and the velocity profiles within the jet. In 
general, flows of this type depend on viscous, surface tension, inertial, and body forces. 
Physical intuition tells us that if the velocity is large and the diameter of the plate is 
sufficiently small, there will be regions wherein viscous forces are not of prime im- 
portance in determining the resulting flow behavior, particularly the free-surface shape. 
The viscous forces, in this case, will be confined to a thin boundary layer on the plate 
that originates from the stagnation point. The location of the stagnation point, O, is 
shown in figure 1. The jet or nozzle radius is Rq, and the distance between the plate 
and the nozzle is H. A cylindrical coordinate system (r,z) emanating from the stagna- 
tion point is chosen. An order-of- magnitude analysis will permit the governing equa- 
tions to be simplified so that an analytical solution can be attempted. For axisymmetric, 
isothermal, incompressible steady flow under weightless conditions, the governing 
equations in cylindrical coordinates for a Newtonian fluid can be written as 

Continuity: 


I JL ( ru ) + = 0 

r 9r 9z 


Momentum: 

r Component, 


(vL — + V — \ = - ^2 + JLL J 

fa 

(ru) 1 

♦ £s] 

V 9r 9z/ 9r 

l 3r 

r 9r 

dz^J 


( 1 ) 


( 2 ) 
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z Component, 


P 


( u 


Jv + v Jv\ = _ JE + 

3r 3z/ 3z 



(3) 


(All symbols are defined in appendix F. ) 

Boundary conditions are required on the flat plate, along the axis of symmetry, at 
the nozzle exit, on the free surface, and at r = L. 


On plate: 


Along axis of symmetry: 


u = O'! 


v = OJ 


■ on z =0, all r 


u = 0 


iz = o 

3r 


> on r =0, all z 


(4) 


(5) 


At nozzle: 


V = -tO 

u = 0 


■on 0 < r < Rq, z = H 


(6) 


At r = L: 


u =u(zf| 

Von r = L, 0 < z < f(L) (7) 

v = v(z)J 

On the free surface, denoted by z_ = f(r), two boundary conditions are required 
since the free- surface position is an unknown to be determined as part of the solution. 
The details of the calculation for the boundary conditions along the free surface are 
given in appendix A (see eqs. (A10) and (A14)). 
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On free surface: 


1/2 2\ a d i 

— (u + v ) - — — | 

2 pr dr 


df 

dr 




■y 


pR 0 


on z = f(r) 


( 8 ) 


and 


-u — + v = 0 on z_ = f(r) (9) 

dr s 

In equation (4), the no- flow and no- slip boundary conditions are applied at the wall. 
Equation (5) is a statement involving the known geometrical symmetry of the problem; 
equation (6) imposes an initially uniform velocity profile on the incoming jet. Equa- 
tion (7) simply states the velocity distribution as the liquid leaves the control volume. 
Equation (8) is a statement of conservation of mechanical energy along a streamline; 
equation (9) states that the normal velocity component on a streamline is zero. The 
second terms on the left and right sides of equation (8) are the contribution of surface 
tension to the mechanical energy balance. 

The solution of the problem can be greatly facilitated by simplifying equations (1) 
to (8). Specifically, the method of obtaining the minimum parametric representation of 
a problem will be used to simplify the governing equations. This method is described 
in detail by Krantz (ref. 38) and is the most systematic approach for scaling the govern- 
ing equations. The initial step in the minimum parametric representation method is to 
form dimensionless variables by introducing characteristic scale factors for all de- 
pendent and independent variables. The unknown scale factors are defined as Uq, Vq, 

i> z 0 ’ p 0> and f 0* 

Dimensionless variables are now defined as 

u* = JL, v * =— , r* , z* =— , p* = -2_, f* =— (10) 

U 0 V 0 r 0 z 0 p 0 f 0 

Introducing these dimensionless variables into the differential equations and boundary 
conditions and arbitrarily making the coefficient of one term in each differential equa- 
tion and boundary condition equal to 1 result in 
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Continuity: 


0v* [ U Q Z Q 1 d 
dz* VqTq r* 0r* 


(r*u* ) = 0 


(ID 


Momentum: 

r Component, 


pU 0 z 0 u * _9u^_ + pV 0 z 0 y * 3u* 
TqP 9r* M dz* 


2 2 

P0 Z 0 gp* [ z 0 d 

r 0 pU 0 ar * r| 9r * 


l d 
r* gr* 


(r* u* ) 


[ 9 2 u* 
dz* 2 


(12) 


z Component, 


U 0 Z 0 dy* , T . SV* p 0 3P* | M2 0 1 S_/ r » (13) 

r 0 v 0 3r * az * pV 2az * r 2 pV 0 r * ar *' ar *' V v 0 3z* 2 

Boundary conditions: 

On wall, 



on z* = 0, all r* 


(14) 


Along axis of symmetry, 


u* = 0) 


lonr* = 0, all z* 

gr* 


(15) 
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At nozzle, 


* 'V 

v* = — — 


J>0 s r* 


u* = 0 


J 




(16) 


At r* = L/r 0 , 


u* = u* (z* f 
v* = v* (z* )_, 


•on r* = L/Tq, 0 < z* s f(L)/z^ 


(17) 


On free surface, 


i /u. 2 ^ + V* 21 

V 2 





2 A p r o v o 


f*f r 


on z* =- 


(18) 


and 


,* _ u * df*^ 


Vo 


dr* 


= 0 


on z* - 
s 


f*f r 


(19) 


The scale factors must now be determined. This is done by setting some of the re- 
sulting dimensionless groups in the equations and boundary conditions equal to zero or 
1. The groups chosen depend on the physical conditions for which the equations are 
being scaled. In what follows, the formalism of Krantz (ref. 38) is strictly followed. 
Characteristic lengths are usually determined from the dimensionless groups generated 
by the boundary conditions. Characteristic times, velocities, etc. , are determined 
from dimensionless groups generated by the differential equations. The guidelines in 
determining the unknown scale factors are 

(1) Do not introduce any mathematical contradictions. 

(2) Do not violate physical intuition. 
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BOUNDARY CONDITIONS 


If we examine the boundary conditions, it is apparent that the following dimension- 
less groups are introduced: 


2 

Y R 0 u L f(L) U 0 



of. 


2 2 


PVfr 


cro 



a f _0 } JVo 

PRqVq z o v o r o 


It is known that v has the range 0 to ~Y, r has the range 0 to L, z has the range 
0 to H, and f(r) also has the range 0 to H. Therefore, setting 



_L 

r 0 



( 20 ) 


implies that 


A 



r 0 



( 21 ) 


and yields two of the six unknown scale factors. Then v* and r* are said to be of the 
same or smaller order of 1. Similar statements can be applied to the remaining vari- 
ables. This ordering procedure is discussed thoroughly by Lu (ref. 39). 

Some of the remaining dimensionless groups cannot be set equal to 1 or zero without 
introducing contradictions. Setting Vq = •f' and = L into these ratios and assuming 
that 

f Q = z Q (22) 

gives the following meaningful ratios: 

_H ^ 02 0 Z Q U Q Z Q 

z 0 "V 2 PT 2 L 2 L 2 
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Setting the second or fourth ratio equal to 1 or zero would violate physical intuition. 
Therefore, the ratios to be considered are 

gz o u o z o 

pT 2 L 2 ’ 

At this point an attempt was made to set H/z^ equal to 1 such that Zq would equal 
H. This seemed logical because 0 to H was the range of z. However, this leads to 
some confusing results in terms of the physics. For a given flow condition, it is argued 
that for a certain (minimum) value of H up to H = 00 , the flow pattern in the vicinity of 
the plate is not expected to change. This is shown schematically in figure 2. This ar- 
gument has been experimentally verified and is discussed at some length in the section 
EXPERI MENTATION . The major point is that H cannot be a characteristic length in 
the problem either with reference to Zq or to fQ. This leaves two remaining ratios 
from consideration of the boundary conditions: ctZq/p^^L^ and UqZq/V'Tj. Accordingly, 
all possible information from the boundary conditions has been obtained. Two of the six 
scale factors and a relation between two others have been determined. The governing 
equations must now be examined. 



GOVERNING EQUATIONS 


Continuity Equation 

From the physics of the problem, it is known that mass must be conserved. If the 
dimensionless derivatives 3v*/3z* and 

— — (r* u* ) 

r* 3r* 


are to be of the same order of magnitude, it is required that 



Solving for U Q gives 


U Q =^ (25) 

z 0 

This, of course, is an equation relating two unknowns, Uq and Zq. The same informa- 
tion could have been obtained by setting the second of the two ratios remaining from the 
consideration of the boundary conditions equal to 1. 


Momentum Equations 

It is the objective of this analysis to define that portion of the flow for which an in- 
viscid solution is valid. For this case, the pressure forces are balanced by the inertial 
forces. This fact allows us to determine the scale factor for the pressure Pq. If the 
dimensionless pressure gradient in equation (13) is to be the same order of magnitude as 
the dimensionless inertia term, 


pv l 


(26) 


it implies that 


p 0 = pV 2 = pf 2 


(27) 


Thus, the characteristic pressure is the stagnation pressure. 


Physics of Zq 

The remaining unknown to be determined is Zq. At this point some physical argu- 
ments are necessary. Recall that it has been shown that H cannot be considered as a 
characteristic scale factor for the problem at hand for reasonably large values of H. 
However, the free-surface shape is expected to change as Rq varies (fig. 3). There- 
fore, from physical considerations, this suggests that characteristic values of Zq vary 
as Rq. By defining 


Zq - Rq 


(28) 
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Uq can be found from equation (25): 


U 0 = 


•f"L 


(29) 


In summary, the following has been determined: 


v 0 = ^ 


r 0 = L 


z o ~ R 0 


f 0 - R o 


u o = 


Po 


£L 

R o 

= pr 2 


These scale factors can now be substituted into the governing equations to obtain the 
minimum parametric representation of the problem. The results are as follows: 


Reu* — — + Re v* Su * 


9r* 


0Z* 


-Re»iPl + pT-^ 

V L / Sr* \ L / Sr* 


— — — (r* u* ) 
r* Sr* 


S 2 u* 

Sz* 2 


(30) 


Re u* 3vl + 
Sr* 


R ev*iZi = -Re-§El + ^ -i-Mr* i^lV 

Sz* Sz* \L / r* Sr* V Sr*/ 


S 2 v* 

Sz* 2 


(31) 


where Re is the Reynolds number, defined as 


Re = 


P^R, 


M 


(32) 
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The equations are now in the form from which it can be determined what the conditions 
must be such that viscous forces are not significant. The major parameters in this 
problem are Rq/L and Re - Since all the starred terms in equations (30) and (31) are 
of order 1 or smaller, only the coefficients of the individual terms need be considered 

O 

in simplifying. Considering equation (30), it can be seen that since (Rq/L) « 1, both 
the inertial and pressure forces -will be an order of magnitude greater than the viscous 
forces, provided that 


Re >> 1 
Re(R 0 /L) 2 >> 

2 

From equation (31), since (Rq/L) « 1, no new information is obtained. The govern- 
ing equations will be reduced to a simplified form of Euler's equations of motion 


) v * 9u* _ _ 

/ R 0\ 2 9p* 

(34) 

9r* 

9z* 

yL/ 9r* 


* 3v* * 3v* 

u * + v * 

9p* 

(35) 

dr* 

9z* 

9z* 



RESULTS 


If the viscous forces are restricted to be at least two orders of magnitude smaller 
than the inertial or pressure forces, there results 


Re > 100 1 


Re 



> 

> 100 

J 


(36) 


Since Rq/L is less than 1 by definition, the coefficient to consider is the second one in 
equation (36), since this will be the limiting one. Under the restrictions Re ^ 100 and 
Rq/L < 1, the Euler's equations of motion are obtained. The equation Re (Rq/L) = 
100 is shown graphically in figure 4. The curve in figure 4 separates the inviscid re- 
gion from the viscous region. Physical understanding of the problem is made clearer 
by referring to this plot. The higher the incoming jet Reynolds number becomes, the 
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thinner will be the boundary layer at some fixed radial position from the stagnation 
point. As shown in figure 4, at lower values of the ratio Rq/L (perhaps obtained by 
increasing the disk radius L), higher Reynolds numbers are required in order to avoid 
viscous influence. Finally, within the viscous region, the boundary layer will grow to 
a radial length that is equal to the free-surface height. The effective design of an ex- 
periment is now possible so that the flow can be considered essentially inviscid. 


EXPER IMENTATION 
APPARATUS AND PROCEDURE 
Test Facility 

The experimental investigation was conducted in the Lewis 2. 2-second drop tower. 
The exact specifications of the facility, the mode of operation, and the release and re- 
covery systems are described in detail in appendix B. The drop tower provides a 2.2- 
second weightless environment in which to conduct tests. 


Experiment Package 

The experiment package used to obtain the data for this study is shown in figure 5 . 

It consists of an aluminum frame in which were mounted the jet reservoir, the disk, the 
16- millimeter high-speed motion picture camera, the liquid supply tank, the backlight- 
ing scheme, and the batteries. The major functions were controlled by onboard se- 
quence timers. 

A diagram indicating the manner in which the flow system operates is shown in 
figure 6. This was a pressure- controlled system in which the flow was begun by open- 
ing the solenoid valve. Prior to the drop, liquid was contained in the line between the 
liquid supply tank and the jet reservoir. In addition, the jet reservoir was completely 
filled with the test liquid. 

A schematic drawing of the jet reservoir, which was fabricated of acrylic plastic, 
is shown in figure 7. The critical feature of the jet reservoir is the 45° tapered section 
that leads to a circular nozzle of diameter dg. The range of nozzle radii studied was 
from 0.25 to 0. 75 centimeter. In addition, the transition between the circular nozzle 
and the tapered section was rounded smooth. This construction prevents boundary layer 
buildup and allows the liquid jet to exit from the reservoir with a nearly uniform velocity 
profile. Based on the experimental results of reference 40, this is more than sufficient 
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to ensure a uniform exiting velocity profile over the range of Reynolds number studied. 

Sharp-edge disks, also fabricated of acrylic plastic, were mounted above the jet 
reservoir by means of a threaded rod such that the flat surface of the disk was at right 
angles to the impinging liquid jet. The radii L of the disks were 1. 0 and 1. 5 centi- 
meters. A schematic drawing of the disks is shown in figure 8. 


Test Liquids 

Two test liquids were used, anhydrous ethanol and trie hlorotrifluoroe thane (Freon 
TF). Their properties at 20° C are listed in table I. No attempt was made to correct 
the fluid properties for temperature changes. Both test fluids possess a nearly 0° static 
contact angle on an acrylic plastic surface. However, this was not the reason why they 
were chosen as test fluids. They were chosen because of their relatively low viscosity 
and their availability. 


Test Procedure 

Prior to a test run, the jet reservoir, the disk, and the supply tanks were cleaned 
ultras onically with a mild detergent. After these parts were rinsed with methanol, they 
were dried in a warm-air dryer. The supply tank was subsequently filled with the test 
liquid, and the jet reservoir was filled by pressurizing the supply tank. This procedure 
eliminated air bubbles from the lines and ensured accurate flow rates. After the jet 
reservoir was completely full, the supply tank was sealed and two accumulator bottles 
(not shown in fig. 5) were pressurized with gaseous nitrogen to a predetermined value. 
The accumulator bottles were designed to be of such a volume that flow rate and pres- 
sure remained essentially constant during a test. Thus, jet velocities were determined 
by measuring the change in liquid level in the supply tanks during a normal- gravity 
calibration run. It was assumed that the flow rate obtained in this fashion would remain 
the same in weightlessness. 

Electrical timers on the experiment package were set to control the initiation and 
duration of all functions programmed during the drop. The experiment package was 
then balanced and positioned within the prebalanced drag shield. The wire support was 
attached to the experiment package through an access hole in the shield (fig. 29(a)). 
Properly sized, spiked tips were installed on the drag shield. Then the drag shield, 
with the experiment package inside, was hoisted to the predrop position at the top of the 
facility (fig. 28) and connected to an external electrical power source. The wire sup- 
port was attached to the release system, and the entire assembly was suspended from 
the wire. After final electrical checks were made and the experiment package was 
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switched to internal power, the system was released. After the test was completed, the 
experiment package and drag shield were returned to the preparation area. 


EXPERIMENTAL RESULTS 
Effect of Nozzle Height 

It was desired to conduct the zero-gravity tests in regions where the nondimensional 
distance between the nozzle and the disk H/R q was not of importance to the problem. 
Several tests were conducted that showed no effect on the free- surface shapes if the ratio 
of nozzle height to jet radius was greater than 6. Therefore, all zero-gravity experi- 
ments were conducted under this constraint in order to eliminate H/Rq as a variable. 
This experimentally determined fact supports the argument made in the section ORDER- 
OF-MAGNITUDE ANALYSIS that H could not be a scale factor for the axial coordinate. 


Steady-State Flow Patterns 

The approximate steady-state flows in weightlessness for three different jet veloci- 
ties are shown photographically in figure 9. The direction of flow of the liquid jet is 
vertically upward. The threaded rod and bolt shown in the figure make up the disk holder 
that connects the sharp-edge disk assembly to the rig frame. The jet velocity increases 
from left to right in the figure. Three distinct classifications of flow pattern were ob- 
served to occur and are labeled in figure 9 as surface tension flow, transition flow, and 
inertial flow. Surface tension flow (fig. 9(a)) is defined as that flow in which the liquid 
flows completely around the disk, with no separation at the disk edge. In transition flow 
(fig. 9(b)), surface tension and inertial forces are both important. Transition flow is 
defined as flow in which the liquid separates from the disk and the resulting liquid sheet 
either collects upon itself and forms an envelope or has the tendency to do so. Inertial 
flow (fig. 9(c)) is defined as that flow in which the liquid separates from the disk with no 
liquid returning toward the jet centerline or attempting to form an envelope. 

The flow pattern shown for transition flow is not the flow pattern one would expect if 
steady-state conditions could have been reached. The reason for this is that some liquid 
always travels toward the disk from the point at which the envelope meets. This liquid 
flow strikes the back of the baffle and subsequently disrupts the initial configuration of 
the envelope. The recirculating flow is a strong function of the geometry of the disk 
holder. The surface tension flow is generally slow and, as a result, does not quite 
reach a steady- state configuration on the back of the disk. The inertial flow tests 
(represented by fig. 9(c)) always reached steady state. 
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The experimental tests were conducted for conditions that are in the inviscid region 
of figure 4. (Two exceptions are noted in table n. ) Depending on the particular value 
of Rq/L, there exists a minimum Reynolds number below which the runs can no longer 
be considered as viscous free. The experimental results are given in table II, including 
all the important parameters as well as the flow classifications. One additional param- 
eter, the Weber number, is given in table I. As will be shown in the next section, when 
the Reynolds number is no longer a parameter of controlling interest, the Weber number 

O 

and the ratio Rq/L become important. The Weber number pV Rq/ct, which naturally 
arises from consideration of the boundary conditions, is the ratio of inertial to surface 
tension forces. In the flow category column of table n, S indicates surface tension 
flow, T transition flow, and I inertial flow. Finally, note that the designated flow 
classification for some cases, particularly those bordering transition or inertial flow, 
could easily fit into either category. 


Zero-Gravity Results 

The data contained in table n are shown graphically in figure 10. The curves in 
the figure were faired by hand and separate the various flow classifications. At any 
particular value of the ratio Rq/L in figure 10, the flow classification depends only on 
the system Weber number. Two critical Weber numbers occur at a constant value of 
Rq/L. The lowest critical Weber number separates the surface tension flow from the 
transition flow; the higher critical Weber number separates the transition flow from the 
inertial flow. In addition, the critical Weber number between regimes decreased as 
Rq/L was increased over the range of parameters studied. 


POTENTIAL FORMULATION 

GOVERNING EQUATIONS AND BOUNDARY CONDITIONS IN 
PHYSICAL PLANE INCLUDING SURFACE TENSION 

In the section ORDER- OF- MAGNITUDE ANALYSIS , it was shown that at any par- 
ticular value of Rq/L, if the Reynolds number is sufficiently large, the flow can be 
considered as viscous free. It is further assumed that the jet will continue to remain 
viscous free after leaving the disk. There will be no shear stress between the exiting 
radial jet and the ambient air surrounding it. 
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General Formulation 


Consider the flow of a circular liquid jet impinging normally on a circular disk, as 
shown in figure 11(a), where Rq is the jet radius, L is the disk radius, and H is the 
distance between the jet reservoir and the disk. The incoming jet velocity is given as 
"V, and the initial velocity profile is assumed to be uniform. There will be two free 
surfaces involved. The upper free surface is defined as z g = f^(r), and the lower free 
surface as z g = f 2 (r ) . In addition, a third surface is required for the complete formu- 
lation of the boundary value problem. Initially, this surface was chosen as the straight 
line F^E (fig. 11(a)). However, this proved to be inconvenient since JFE is arbitrary 
and, thus, has no known boundary condition. Instead, z g = fg(r) was chosen, a surface 
along which the velocity potential is constant. Various points in the physical plane have 
been designated with letters ranging from A to G for convenience. In cylindrical 
coordinates, the dimensional governing equations and boundary conditions in terms of 
the primary variables (u,v) are given as follows for a weightless environment and for 
inviscid, incompressible, isothermal, steady flow as 

Continuity: 


A JL (ru ) + ii = 0 

r 3r 3z 

Momentum: 

r Component, 

p ( u Ju + v iE\ = - JE 

\ 3r 3z / 3r 

z Component, 

p ( uix + vixWie 
\ 3r 3z/ 3z 

The following boundary conditions are applied: 

On DC: 


v = 0 on z = 0, 0 2 r s L 


(38) 


(39) 


(40) 


(41) 
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On AB: 


u = 0 0 s r s Rq, z = H 


(42) 


On BC: 


-^ = 0 r=0, 0<z<H 

ar 


On AG: 


A (u 2 + v 2 


)-JLAl 

pr dr 


rfl 






pR ( 


on z g = fj(r) 


(43) 


(44) 


and 


-ufj + v = 0 on z g = fj(r) 


(45) 


On DE: 


- (u 2 + v 2 



= A^ 2 -_2_ onz =f 2 (r) 
2 PR 0 s 2 


(46) 


and 


-uf£ + v = 0 on Zg = fg(r) 


(47) 


On GE: 


u + vf£ = 0 on z g = f 3 (r) (48) 

where fg(r) for now will be assumed to be some specified surface. 

The derivation for the boundary condition along GE is as follows: The velocity 
potential <p is constant. By definition, the velocity vector V is normal to an equipo- 
tential surface. This means that 

V x n = 0 on z g = f g (r) (49) 
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The unit normal to fg(r) is 


n 


f 3* + * 


(50) 


where i is the unit vector in the radial direction, j is the unit vector in the axial direc- 
tion, and 


V = ui + vf 


(51) 


Applying equation (49) results in 


u + vfg(r) 


= 0 


(52) 


From which equation (48 ) follows . 

Introduction of Stokes stream function . - The primary variables contained in the 
governing equations are the scalar velocity components u and v. The fact that two 
functions are required to describe one vector field is cumbersome. As shown in the 
theory of hydrodynamics, the number of functions can be reduced for several important 
cases, one of these being axisymmetric flow. A function \p , defined as the Stoke's 
stream function, can be introduced, which automatically satisfies the continuity equa- 
tion. With the additional requirement of irrotationality, the governing equation in terms 
of the Stokes’s stream function will result. 

According to Chan (ref. 27), the Stokes's stream function is a mathematical device 
used to describe the flow. It has the following properties: First, when the stream 
function is set equal to a constant, it results in different annular stream surfaces in 
axisymmetric flow. Second, when it is differentiated properly, it yields the velocity 
components. Third, the difference between the values at two adjacent stream surfaces 
is related to the flow rate between those surfaces. 

Starting with equation (38), the continuity equation for axisymmetric flow, the fol- 
lowing substitutions are made: 


u = -i^ 

r 3z 


(53) 
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and 


y =1M. (54) 

r 3r 

Substituting equations (53) and (54) into equation (38) shows that continuity is identically 
satisfied. In addition, the assumption is made that the flow is also irrotational. As a 
result 


curl V = 0 


For axisymmetric flow, this can be written 


3u _ 3v _ o 
3z 3r 


(55) 


(56) 


Replacing u and v in equation (56) by their relations to \f/ results in the governing 
equation for axisymmetric flow in terms of the Stokes’s stream function 




3z‘ 


r 3r 


3r 


(57) 


Thus, one equation, equation (57 ), takes the place of the continuity equation and the two 
momentum equations . As mentioned, continuity is identically satisfied by substitution 
of equations (53) and (54). In addition, substituting equations (53,) and (54) into the 
momentum equations and subsequently eliminating the pressure terms by cross differ- 
entiation show that it too is identically satisfied provided that curl V = 0, which in 
essence is equation (57). 

Derivation of boundary conditions in terms of \j / . - To derive the boundary condi- 
tions in terms of the Stokes stream function \}/ , we need only refer to the physical- 
plane conditions (eqs. (41) to (48)) and the two definitions (eqs. (53) and (54)). 

On DC: 


v = 0 onz = 0, OsrsL 


(58) 


Using equation (54), 
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( 59 ) 


A Ilk = 0 on z = 0, 0 s r < L 
r 3r 


Therefore, 


— = 0 on z = 0, 0 < r =£ L 
9r 


Hence, 


xp # xp(r) on z = 0, 0 s r s L 


and at most. 


ip = xp(z) on z = 0, 0 ^ r < L 


But, z = 0. Therefore, 


(60) 


(61) 


(62) 


xp = Constant = C 

and DC is a streamline, as expected. Arbitrarily choosing the value of this constant 
equal to zero gives 

On DC: 


xp — 0 onz = 0, OsrsL (63) 

Since DC, BC , and DE are all part of the same streamline, they must all have the 
same Stoke's stream function (xp = 0). As indicated in equation (63), this value is ar- 
bitrarily chosen as zero. 

Now, examine the surface AB . 

On AB: 


u = 0 0 < r s R q , z = H (64) 


By using equation (53 ), it can be seen that 

u=-AAAl = 0 0<r<R 0 , z=H (65) 

r az 
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and, therefore, 


-^ = 0 0 < r ^ R 0 , z = H 


dz 


( 66 ) 


and generally 


xp = xp(r) 0 < r < R Q , z = H 


(67) 


However, we know on AB 


v = -Y 0 s r < R n , z = H 


(68) 


so that 


■f =-^ 0 < r < R ft , z=H 

r 3r u 


(69) 


Integration yields 


xp = - + C 0 < r < R n , z = H 

2 u 


(70) 


The value of C can now be determined by evaluating this last expression at r = 0. 


* 


r=0 


= C 


(71) 


But along BC, xp - 0. Therefore C = 0 at r = 0 and 


xp = - 0 < r s R n , z = H 


(72) 


On AG: Since AG is a line of constant xp , applying equation (72 ) yields 


xp =- 


YR ( 


(73) 


In addition, equation (44) applies. Substituting equations (53) and (54) into equation (44) 
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yields 




On GE: Equation (48) become 

( r ) = 0 on z = fq(r) (76) 

dz dr 6 s 6 

The second boundary conditions on the free surfaces, namely, -uf^ + v = 0 and 

-ufA + v = 0, are fully equivalent to the specification of the value of the Stokes’s stream 

" 2 
function along those surfaces, ^ = 0 on DE and \J/ = - t^Rq/2 on AG. Since i// = 

xp{r,z), it can be expanded and d\p is equal to zero along AG and ED . 

Nondimensionalization of governing equations and boundary condi tio ns i n physical 

plane . - The governing equation in terms of the Stokes’s stream function (eq. (57)) and 

the boundary conditions are now put into dimensionless form by introducing scale fac- 

2 

tors. Let the scale factor for the stream function be -^Rq and Rq be the scale 
factor for lengths. Let dimensionless quantities be represented by stars, (i.e. , if/* is 
dimensionless). The results of this manipulation are as follows: The governing equa- 
tion is 



- A. Ml. = 0 

r* 3r* 


(77) 


The boundary conditions are 
On DC : 

xf/* = 0 on z* = 0 0 < r* s L/R q (78) 
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On AB: 


On BC: 


On AG: 


and 


p* = - r* 2 0 < r* < 1, z* = H/R n 

2 


<//* = 0 r* = 0, 0 < z* < H/R Q 


\I/ * = — on z* = f f (r* ) 
2 si 


(79) 


(80) 


(81) 



Three parameters appear in the specification of boundary conditions. They include 
the Weber number We and the dimensionless length ratio L/Rq. Recalling the argu- 
ments in the ORDER- OF-MAGNITUDE ANALYSIS , section H/Rq is not really a param- 
eter if it is larger than some minimum value. The dimensionless velocity components 
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can be calculated from the following expressions: 


u 

i dp* 

(86) 

Y " 

r* dz* 

v _ 

1 3^* 

(87) 

Y 

r* 9r* 


The procedure for numerical solution of the boundary value problem using finite 
difference methods as set up in dimensionless form in the physical plane would be as 
follows: Initially, realistic variations for f|(r), f|(r), and f| (r) must be assumed. 
Using one of the two given boundary conditions on AG and DE(i.e. , p* =1/2 and 
ip* = 0), solve for ip * . With the initial solution for p* , check the validity of the 
second of the two boundary conditions on AG and GE . If the boundary conditions are 
not satisfied, new variations in f| , f^ , and f| must be assumed. A second iteration 
to ip* must be obtained and so on. One of the serious drawbacks of this outlined iter- 
ation scheme is the lack of knowledge concerning how to update assumed values of f| , 
f|, and fjjj based on previous solutions. In other words, there is no logical way in 
which to make changes to the shape of the initially assumed control volume. 


Surface- Tension- Dominated Model 

For Weber numbers between 5 and 30 (depending on the ratio Rq/L), experimental 
data (fig. 10) show that the resulting steady- state flow pattern is surface tension domi- 
nated (fig. 9(a)). By previous definition, in surface tension flow the liquid flows com- 
pletely around the disk with no separation from the edge. It is the intent to model this 
flow in order to solve for the theoretical free-surface shapes and velocity profiles. 
Assuming axisymmetry, the physical plane model is shown in figure 11(b). In the model, 
at some cross section far downstream, the flow is assumed to approach the initially 
uniform flow it possessed at AB. The exiting plane is denoted by GE in the model. 
Both C and C* are located at r = 0, z = 0. (The plate is assumed to be infinitely 
thin. ) However, C is located on top of the plate, while C’ is on the bottom. The free 
surface fj(r) is not assumed to possess mirror symmetry about the z = 0 position. 

In cylindrical coordinates, the governing equations and the boundary conditions in terms 
of the primary variables (u, v) are given as 


29 



Continuity: 


I-L (ru ) + 8 t . 0 

r 3r 9z 

Momentum: 

r Component, 

p fuiH + v iH\ = -lP 

\ 9r 9z/ dr 

z Component, 

p/uix + vixU-iE 

V 9r 9z/ 3z 


The following boundary conditions are applied: 


On DC: 


On DC’: 


v - 0 on z = 0, O^r^L 


v = 0 on z = 0, Osr <L 


On AB: 


u=0 0sr<R 0 , z = H 


On BC: 


On C'E: 


— = 0 on r = 0, 0 < z < H 
3r 


— = 0 on r = 0, -H s z < 0 
9r 


( 88 ) 


(89) 


(90) 


(91) 


(92) 


(93) 


(94) 


( 95 ) 
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On GE: 


On AG; 


u = 0 0 < r < R Q , z = -H 


1/2 2\ ad 

- (u + v ) — 

2 pr dr 



a 

P R 0 


°n z g = fj(r) 


(96) 


(97) 


and 


-uf^ + v = 0 on z g = f j(r) 


(98) 


Equations (97) and (98) represent two distinct pieces of information concerning z = f^(r). 

Direct substitution of equations (53) and (54), the relations between the velocity 
components and the Stoke's stream function, into the governing equations and the bound- 
ary conditions for the surface tension model results in the desired formulation after 
nondimensionalization. Similar to the general formulation in the last section, the vari- 
ous lengths in the problem are scaled with respect to Rq and the scale factor for the 
Stoke's stream function is -VRq. The equations are 


5 2 ip* ] 9 2 ^/ * _ 1 9;f/* _ o 

9z * 2 9r* 2 r * 9r * " 


(99) 


The boundary conditions are 


On DC: 


On DC': 


On BC: 


\j/* = 0 on z* = 0, 0 < r* < L/R q 


ty* = 0 on z* = 0, 0 < r* < L/R q 


\f/ * = 0 on r* = 0, 0 < z* s H/Rq 


( 100 ) 


( 101 ) 


( 102 ) 
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On C'E: 


xf/* = 0 on r* = 0, -H/Rq < z* < 0 


On AB: 


xfz * 


1 *2 
= — r 

2 


on 0 s r* < 1, z* = H/Rq 


On GE : 


On AG: 


xp* = -r* 2 on 0 < r* < 1, z* r -H/R n 

2 U 


xp* =- on z* = f*(r*) 
2 s 


(103) 


(104) 


(105) 


(106) 


and 



on z* = f*(r*) 


(107) 


INVERSE -PLANE FORMULATIONS 

The procedure for solving for the free-surface shapes in the physical plane has 
already been outlined. The difficulties encountered when making adjustments to the free 
surfaces between iterations and the lack of a logical manner in which to make the ad- 
justments have been cited. This would be a time-consuming task even in the absence of 
surface tension. A computerized scheme is sought that offers the possibility of achiev- 
ing the free-surface results with a minimum of computer iteration time and user inter- 
action. 
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Transformation Formulas 


An alternative approach to the physical- plane solution is discussed in detail by 
Jeppson (ref. 14). In his article, Jeppson discusses a transformation technique into 
what is termed the inverse plane. The coordinate axes in the inverse plane are the 
velocity potential q> and the Stoke's stream function . The advantage of using the 
cp\f/ space, or inverse plane, is that the free surface lies along a line of constant xf/ 
and is, therefore, at a known position. The governing equation, however, is no longer 
linear as it was in the physical plane. 

Introduction of v elocity potential . - A scalar function defined as the velocity poten- 
tial has been used previously in the specification of the boundary GE (the exiting plane). 
The continuity equation for steady incompressible axisymmetric flow is given by equa- 
tion (39). The flow is also as stun ed to be irrotational such that the condition given by 
equation (56) is also valid. Equation (56) implies that there exists a scalar potential 
function <p such that 


V = grad cp 


from which it follows that 


u=i£ 

9r 


and 


v 


9z 


(108) 


(109) 


( 110 ) 


Substituting equations (109) and (110) into equation (38) yields the governing equation for 
steady axisymmetric flow in terms of the velocity potential 


9 2 <? + 1 9£ 
9z 2 r 9r 



(HD 


Relationships bet ween ijy and cp . - It can now be stated that the relations between 
the velocity potential and the Stokes’s stream function are 

U=^- = - (112) 

9r r 9z 
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and 


v = 12- riil 
8z r 3r 


(113) 


Inverse-plane formulation . - To formulate the problem in the inverse plane, r and 
z must be known as functions of xp and <p. In other words, by working in the inverse 
plane we are attempting to reverse the roles of the dependent and independent variables. 
The required relations are obtained by noting that if xp = xp( r, z) and <p = <p( r,z), there 
then exist the inverse functions r = r(<p , xp ) and z = z (cp ,xp): 


and 


r 2 8 2 z + 8 2 z _ _2 _8z_ _8z_ 


dip* d<p‘ 


dip dcp 


(114) 


8 2 r + J_ j?r_ _ J_ / 9r\ 2 + 1 /9r\ 2 _ Q 
dxp 2 r 2 dcp 2 ' A**' r W/ 


(115) 


Details of the derivation of these equations are contained in appendix C and were obtained 
from reference 14. 

For further discussion, equation (114) is referred to as the z equation and equa- 
tion (115) is called the r equation. Both these equations are nonlinear, with the non- 
linearity in the z equation also involving r. Also, both these equations are of the 
elliptic type, as shown in reference 14. This means that boundary conditions are re- 
quired for all boundaries in the flow region. A solution to the problem begins with a 
solution to the r equation (eq. (115)) since it contains only one unknown, r. 

One final point involves the fact that the variables discussed in the inverse trans- 
formations are all dimensional. Nondimensionalization of these equations similar to 
the nondimensionalization done in the physical plane allows the same equation to be re- 
covered. Let 
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z 

z* = — where Zg = Rg 
z 0 

r* = — where r 0 = Rg 
r 0 

tp* = JL. where xp n = -^R« 
*0 

where <Pq = -^Rg 

^0 J 


(116) 


These substitutions yield exactly the same governing equations and transformation 
formulas in starred notation. In other words, 


r * 2 o z* t 


0 2 z* 


dif/ 


* 2 




*2 


= -2 


9z* az* 
dxp* 3<p* 


(117) 


and 


+ JL/J£Lf =o (ns; 

3^*2 r* 2 dcp* 2 r* 3 ' d< P*' r * W*/ 

Infinite Flat Plate Excluding Surface Tension 

The first problem to be examined in which the inverse transformation is used is the 
flow of a circular liquid jet normal to a flat plate (fig. 12). At this point, the starred 
notation to designate dimensionless quantities is dropped. It is assumed that all quan- 
tities appearing henceforth are dimensionless unless otherwise stated (i.e. , r, z, xp, 

<p are now dimensionless). For the case in which surface tension is excluded, some 
of the physical- plane boundary conditions must necessarily change. The model de- 
scribed in equations (77) to (85) is for the general case of a finite disk in which surface 
tension effects are important. The boundary conditions that remain valid for the 
infinite-plate case excluding surface tension are equations (79) and (80), with equa- 
tion (78) modified to read "on EC . " For the case of an infinite flat plate, since the 
free surface DE of figure 11(a) no longer exists, the designation D does not appear in 
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the model. Furthermore, for the case in which surface tension forces are unimportant, 
the boundary condition along the upper surface (referring to eq. (82)) can be written as 
follows: 

On AG: 



(119) 


This expression is obtained by allowing the Weber number (the ratio of inertial to sur- 
face tension forces) to approach infinity. The case we are examining is one in which 
the surface tension forces become vanishingly small or the inertial forces become large. 
Finally, the boundary GE will be nearly vertical, implying that the velocity there is 
purely radial. This is tantamount to letting f| ' approach infinity in equation (85). 

Derivation of boundary conditions in inv erse plane for r formulation. - The bound- 
ary conditions in the inverse plane for r formulation are as follows: 


On AB (or (p = Constant): 



in the physical plane, eq. (79) 


r 



in the inverse plane 


(120) 


On BC (or p = 0): 


On EC (or i// = 0): 


r = 0 in the inverse plane 


( 121 ) 


v = 0 

By equation (113), this implies that dcp/dz and d\p/<pr equal zero. In the inverse 
plane (by eq. (C4)), 


0r/a^ = 0 (122) 

On GE (or <p = Constant): It is arbitrarily set equal to zero, cp = 0. It will be assumed 
that GE is nearly vertical, implying that the velocity is purely radial. From equa- 
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tion (85), d\p/dr = 0. In the inverse plane, 


dr/di}/ = 0 since J * 0 


(123) 


On AG (or ^ = 1/2): 



in the physical plane (from eq. (119)) 


From equations (112) and (113), 



= 1 


From equations (C2) and (C4), 



2 

Substituting for J from equation (C9) and then using equations (C6) and (C7) yield the 
final form 



(124) 


The results of these derivations are shown schematically in figure 13(a) in the inverse 
plane. The boundary conditions for the z formulation are also shown in figure 13(b). 

Derivation of boundary conditions in inverse plane for z formulation . - The bound- 
ary conditions in the inverse plane for z formulation are as follows: 

On AB (or <p = Constant): In the physical plane, u = 0, which implies that dcp/dr = 0. 
From equation (C2), 9z/3 \p = 0. This implies that z =z(<p) alone and, since AB is a 
line of constant cp , 


z = Constant in the inverse plane 


(125) 
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On BC (or ip = 0): In the physical plane, u = 0, from which it follows, in the inverse 
plane, that 


-25- = 0 (126) 

dip 


On EC (or ip = 0): 


On GE (or <p = 0): 


z - 0 in the inverse plane 


(127) 


v = 0 in the physical plane 

From equation (113), it follows that dip /dr = 0 and from equation (C3) that, in the 
inverse plane, 


-25- = 0 (128) 

d(p 


On AG (or \p = 1/2): 



in the physical plane 


In the inverse plane we can take over the equivalent expression given by equation (124), 



= 1 


Using equations (C6) and (C7) to eliminate the derivative involving r results in 



(129) 


The solution technique for the infinite -plate case is as follows: The r values can 
be completely found once the governing equation (eq. (118)) is solved. The variable z 
appears nowhere in the formulation (i. e. , it does not appear in either the governing 
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equations or the boundary conditions). Once the r values are found, the z governing 
equation (eq. (117)) can be solved. The z solution requires a prior knowledge of the 
values of r both in the interior region of flow and on the free- surface boundary. 


Finite Plate Excluding Surface Tension 

A second problem to be formulated is concerned with the impingement of a circular 
liquid jet normal to a disk of finite width. For the case in which surface tension effects 
are neglected, the formulation depicted in figure 13 will vary only slightly for the finite 
plate. In fact, all the boundary conditions are identical with the exception of the bottom 
free surface (ED ) and the exiting plane (GE ). 

On ED: 



(130) 


There is an additional change in the boundary condition along GE (or cp = 0), the exiting 
plane of the jet. In the physical plane, from equation (85), 


Ok- - f ' = 0 

dz 3r 6 


By using equations (Cl) and (C3), this can be written 


— + — fl =0 

d<p dcp ^ 


Finally, using equation (C7) to eliminate 3z/3 <p yields 


*L + r**-fU=0 
3 (p d}fy 0 


or 


3 <p vty/\dr 



(131) 


(132) 


(133) 


(134) 
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The boundary conditions for the z formulation in the inverse plane can be obtained in 
a similar manner. Continuing with equation (134) and using equations (C6) and (C7) 
give 

r fiSA- (—)(—)= ® (135 

\3 ipj \d<p/\drj 

The formulation of the problem of flow of a circular liquid jet normal to a finite 
disk in the inverse plane is now complete. The r solution for this case, however, can 
no longer be obtained independently of the z solution due to the exiting jet boundary 
condition. Thus, the problem will necessitate a simultaneous solution of two partial 
differential equations of the nonlinear elliptic type. 


Surface-Tension-Dominated Model - Finite Plate 


The final problem to be examined involves the surface- tension- dominated flow de- 
scribed in figure 11(b). The inverse-plane formulation for the r solution is indicated 
in figure 14(a). The boundary conditions are the same as derived previously, with the 
exception of the free- surface boundary and the exiting jet condition. Along AG in the 
physical plane, 



This will now be put into inverse form. In the derivation of equation (124), it was shown 
that the first bracketed term on the left, 



(137) 


In addition, the term within the radical becomes equation (138) by using equations (Cl), 
(C3 ), and (C7). 
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Therefore, the inverse boundary condition along AG (or \p = 1/2) becomes 


(138) 


1 



which can be written 


1 



lid 

r 

_ 1 . 

1 

We r dr 

J 1 p r /^f + i 

1/ r 2 \3r /dil/J 

2 

We 


1 1 d 

r 2 

_ 1 _ 1 

We r dr 

J(dr/d<pf + r 2 
J \3r/ dtp) J 

2 We 


Multiplying through by 2 and rearranging yields 



1 - 


2 2 1 d 

We We r dr 



(139) 


(140) 


(141) 


There has been no z symmetry assumed for the r formulation, only axisymmetry. 
The r formulation does not involve the variable z. Therefore, the r and z solutions 
can be obtained independently. The z formulation is shown in figure 14(b). 


CENTRAL FINITE DIFFERENCE REPRESENTATION - RE CTANGULAR MESH 

The finite difference operators for the nonlinear r and z elliptic partial differen- 
tial equations resulting from the inversion were put into difference form. It was known 
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from experience that considerable flexibility resulted if the difference equations were 
derived by using rectangular mesh. The finite difference expressions along the free 
surfaces (i. e. , AG and ED ) contain a square root, the argument of which becomes 
negative during iterations. The use of rectangular mesh alleviated this difficulty and 
also allowed us to control the size of the flow region. 


GENERAL FORMULATION 
General Considerations 

The partial derivatives appearing in the r and z governing equations are replaced 
by algebraic central finite difference operators. A complete derivation of the various 
operators is contained in reference 41. The notation used in this report is that shown 
in figure 15. If the finite difference analogy of d/dr is 6/6r, the r derivatives can 
be written at point 0 as follows : 


5r r l - r 3 

6 <p 2 A cp 


(142) 


5r _ r 2 - r 4 
5\J/ 2 Ai// 

5 2 r _ r l + r 3 ~ 2r 0 
d<p2 A (p^ 

5 2 r _ r + r 4 ~ 2r 0 

a . 2 . . 2 

dx}/ Ay/ 


(143) 


(144) 


(145) 


Similarly, the central finite difference representations for the 


z derivatives are 


5z Z 1 ~ z 3 

6 (p 2 A (p 


(146) 


5z z 2 ~ z 4 
8 ^/ 2 


(147) 


42 



6 2 z Z 1 + z 3 2z 0 


(148) 


5 2 z _ z 2 + z 4 ~ 2z p 
6 j // 2 A )// 2 


(149) 


Interior Nodal Points 


If equations (142) to (145) are substituted into the r governing equation (eq. (115)) 
with a = Acp/Ap , the finite difference expression for all interior nodal points is 



r 


+ r 4 )+ r 


2 

0 




A( r i + r 3 ) + -l. 


( r i - 


2cr 


8 of 


(150) 


In addition, if equations (146) to (149) are substituted into the z governing equation 
(eq. (114)), the finite difference expression for all interior nodal points is 



+ z 4 ) + 


1 

2(r 2 Q a 2 + 1) 


( z l 


+ z 3 ) + 



( z 2 " z 4^ z l " z 3^ 

(151) 


where, as in the derivation of equation (150), 


_ A cp 


(152) 


FORMULATION EXCLUDING SURFACE TENSION 

The governing equations and boundary conditions were put into difference form using 
central finite differences. Details are contained in appendix D. 
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Infinite Flat Plate 


The finite difference representation for the infinite flat plate is shown in figure 16. 
The algebraic expressions for the boundaries are derived by simultaneous application of 
the governing equation and the boundary conditions at a fixed point. This application 
involves a fictitious point, f, outside the boundary that is subsequently eliminated. 

Point G represents a special point in the finite difference representation for the z 
solution. The governing z equation (eq. (151)) was applied at point G, which resulted 
in two fictitious points. Then the boundary conditions along both AG and GE were 
applied at point G. This allowed the two fictitious points to be eliminated from the 
resulting finite difference expression. Detailed calculations for the boundaries are con- 
tained in appendix D for all models. 


Finite Plate 


The changes in inverse formulation for the finite- plate problem are now presented. 
For this case, the difference operator along GE is more complicated than in the 
infinite -plate case. The following algebraic expression is valid along GE: 


r n - — (r 9 + t a ) + rj 


— - - (r 2 ~ r 4 ) 2 

„2 8 2 4 


2 of 


2r l + ar 0 ( z 2 - z 4 } 


+ -^(z 2 -z 4 )‘ =0 
8 1 4 

(153) 


On DE, the following expression is valid: 


r 


4 

0 




( r l + r 3 ) 



-1 ''I - ^ 

a 


4a" 


= 0 


(154) 


In fact, both G and E represent special points in the r formulation. However, 
one of these, point E, is specified as a known position (r = Constant). At point G the 
equation to be satisfied in the r formulation is 
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: o - r 0 r 


/a ~ r 0< z 0 ' z 4> 


+ r 2 (z Q - z 4 ) 2 - 5- = 0 (155) 


For the inverse z formulation, the following changes occurred from the infinite-plate 
description: 


On GE : 


rp(z 2 + z 4> 2z l - r 0 a ( r 2 - r 4> + (z 2 ' z 4>[ r 0 a(r 2 ~ r 4 ) ] 
2 U + X) + 1) 4 (r 2 + 


and on DE: 


2 HM 


2z 9 - — 
2 r n 


L l / 4a2 - A . 
‘‘o f a‘ 


- J_ ( Zj - z 3 ) 2 + — ^ < Z 1 + z 3 ) 

a 2 J 2(rgOT + 1) 


1^2 _ J_ (Zi 

( r ° a+ i)’ ° 2 


Both points G and E are special points in the formulation. The following equations 


hold there: 


At point G: 


I 



z 0 


0 




- V &2 ■ r 0^ r 0 - r 4> 2 + z 4 


Z 1 - r Q a ( r Q - r 4> 

2 2 . 
r 0 « +1 


a 


2 1 1 

r o + ~; 


a 2 " r 0^ r 0 ■ r 4 > ( 158 ) 


a 


At point E: 


.2 1 


r 0 + 


V a2 - r ?< r 2 - r o > 2 


0 


Z 1 - r Q a ^ r 2 - r 0> 

2 2 i 

r o a +1 


a 


2 1 ’ 

r o + ~ 


a2 " r 0^ r 2 - r 0 )2 < 159 ) 


a 


Detailed calculations for all the additional boundaries encountered for the finite-plate 
case are contained in appendix D. 


FORMULATION INCLUDING SURFACE TENSION 

The finite difference formulation corresponding to the surface tension model 
(fig. 14) is shown in figure 17. Details of these derivations are given in appendix D. 

The r solution can be obtained independently of the z solution (fig. 17(a)). This im- 
plies that once the r solution has been found, there is no iteration required after the z 
equation has been solved. 


DISCUSSION OF NUMERICAL TECH NIQUES 

Since the governing equations for the r and z formulations are nonlinear in the 
inverse plane, in general, the finite difference operations at the interior and boundary 
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points will be nonlinear (eqs. (150) and (151)). For the finite- flat- plate case, as shown 
in figure 16(a), the solution begins with r = r(i// , <p ). With a knowledge of the r solu- 
tion, the z formulation (fig. 16(b)) is solved for z(\[/ , cp). However, in the case where 
the plate is finite, the r formulation also contains z along the exiting plane GE 
(eq. (153)). As a result, the r and z formulations must be solved simultaneously. 
The surface tension model described in figure 17 also allows the r equation to be 
solved independently of the z equation since z appears nowhere in the r formulation. 

In any case, when solving the r equation, the finite difference representation of 
the problem results in N nonlinear algebraic equations in N unknowns . A variety of 
methods were applied to obtain a solution to the simultaneous nonlinear equations . 

These included Lieberstein's extension of Young's work on overrelaxation to nonlinear 
elliptic partial differential equations (ref. 42) and the familiar Newton- Rap hs on method. 
None of the methods were successful in obtaining a convergent solution. 

The technique suggested by Powell (ref. 43) led to the method used in this report 
to obtain solutions . Basically, Powell developed a subroutine that was essentially a 
"compromise between the Newton- Raphson algorithm and the method of steepest de- 
scents. ” Powell describes a FORTRAN subroutine for solving the nonlinear set of 
equations 


f K (Xj,X 2 , . . .,X N ) = 0 K = 1, 2, . . . , N (160) 

The objective is to minimize the function 

N 

F(X 1 ,X 2 ,. . .,X N )= £ [f K (X l5 X 2 ,. . .,X N )] 2 (161) 

K=1 

As with many iteration schemes, initial guesses are required for Xj^. This particular 
algorithm has an advantage in that the initial guessed values do not have to be close to 
the exact solution. The computer program for the r solution contains the main pro- 
gram and three subroutines. The subroutine EQUATIONS is the one supplied by Powell. 
The user supplies the subroutine MATINV, which inverts the matrices, and the subrou- 
tine CALFUN, which contains the nonlinear functions F(Xj^). A knowledge of the r 
solution makes the z formulation explicit in the unknown Zq at each nodal point. As 
a result, a Gauss- Siedel linear iteration scheme was employed to obtain the solution. 

For details of the subroutines, see reference 43. 
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NUMERICAL RESULTS 


FORMULATION EXCLUDING SURFACE TENSION 
Infinite Flat Plate 

Initially, the infinite- flat- plate problem (Weber number = °° ) was numerically 
solved by using a very coarse mesh (fig. 16(a)). The total stream function \}/ was 
divided into five equal parts. Recall that \pQ - 1/2 and = 0, which meant that the 
parameter a = was set equal to 0. 1. In addition, the <p - axis was divided into eight 
equal parts. This division was purely arbitrary. One consideration was that there 
would be at least two vertical lines of constant cp between points C and B. This re- 
sulted in a total of 39 unknown nodal points for the r formulation in figure 16(a) and 
41 unknown nodal points for the z formulation in figure 16(b). The value of D = 1/a , 
where a = A^c/A^/, was set equal to 0. 0204 (a = 7). In figure 16(a), the value of C 
along GE was arbitrarily chosen as 4. 0; in figure 16(b) the value of Kcp for z along 
AB was chosen as 3. 315. The major reason for the assignment of these two values for 
r and z was to be able to compare the results with the semianalytical work for the 
infinite-flat-plate case of Schach (ref. 29). Actually, the major variable choice in the 
entire program is D = 1/a . Basically, D is a measure of how large the flow system 
is since rectangular mesh is being used, that is, a measure of A cp and -j. The 

only concern is that is not too small. In that event, the extent of the solution 

in physical space would be such that ^ total anc * i^ e AB boundary condition would be 
forced too close to the plate to result in a realistic solution in the total flow domain. If 
A <p is chosen too large, that is, a is chosen too large or D is chosen too small, all 
that is lost is accuracy due to larger mesh spacings. In other words, the AB boundary 
is now too far from the plane. For the coarse- mesh, infinite- flat-plate problem, the 
results of the numerical solution are given in appendix E. For the computer listings 
and final output, see reference 44. The solution to the formulation (fig. 16(a)) required 
1098 calls of the subroutine used to solve the simultaneous nonlinear equation. The sum 
of the squares of the error to the solution was reduced to 0. 00357 at the last iteration. 

The method applied to obtain the z solution was a simple linear iteration technique 
since equations explicit in Zq can be derived both in the interior and along the boundary 
points. There were 249 iterations required for the z solution. Changes between the 
248th and 249th iteration occurred in the fifth decimal point. In general, iterations 
throughout the numerical computations were continued until no significant changes in the 
solutions resulted. 

A physical-plane description of the infinite- flat-plate solution is presented in fig- 
ure 18. Only two of the four internal streamlines are shown in the figure. Curves were 
faired through the available calculated nodal points (r, z) by a best-fit process. There 
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appeared, at the onset, the question of whether the coarse mesh employed could suffi- 
ciently describe the flow. When applying boundary conditions at the free surface, a 
larger number of nodal points are desirable. As a result, the existing mesh was 
doubled. The total i f/ was divided into 10 equal parts such that a = A\[/ = 0. 05, and 
the total cp was divided into 16 equal parts. The final values for the coarse-mesh 
solution (r, z) were used to make initial guesses of the values for the fine- mesh solution. 
The r solution required the simultaneous solution of 159 nonlinear algebraic equations. 
The z equation, again being explicit in Zq, resulted in 164 unknown values of z. The 
results from the computer program (results shown in appendix E) are presented in fig- 
ure 19 as computer plot. The computer connected the nodal points with straight line 
segments. In general, there is very little difference between the fine-mesh and coarse- 
mesh solutions. The fine- mesh solution required the extended- storage -space option on 
the 1106 Univac computer. The sum of the squares for the final r solution was reduced 
to 0. 024 after approximately 1000 calls of the subroutine CALFUN used to solve the 
simultaneous nonlinear equations. For this case, 1000 passes through the flow field 
were required to obtain a satisfactory z solution. 


Finite Plate (Inviscid, Excluding Surface Tension) 

As mentioned in the section DISCUSSIO N OF NUMERICAL TECHNIQUES , the finite- 
plate formulation also contains z on the exiting jet surface (eq. (154)). Initially, values 
for z along GE were assumed. These values were used to compute an r solution. 

The computed r solution used 161 calls of the subroutine CALFUN to reduce the sum of 
the squares of the residuals at the nodal points to 0. 00066. With this r solution, the 
linear iteration technique was used to calculate the complete z solution. Since the 
computed z solution resulted in refined approximations to the values of z along GE , 
changes could then be reflected in a new r computation. 

For the finite-disk case in which the ratio of the radius of the liquid jet to the radius 
of the disk was 1/2, a second r calculation did not change when the z values were up- 
dated. A curve was faired through the calculated nodal points (fig. 20). Only two of the 
four available streamlines are shown in the figure. No attempt was made to refine the 
solution by completely doubling the number of vertical and horizontal grid lines. For 
this particular case, since no comparison with any existing analytical techniques existed, 
a more convenient value of 3. 3 was chosen for K cp . The value of D employed in the 
solution was 0. 0138. The only specification along GE was that z at point G be set 
equal to 0.25. The calculated r and z values at each nodal point are given in appen- 
dix E. 

The same method was used to numerically compute the finite-plate case in which 
the ratio of the jet radius to the disk radius was 3/4. There were 38 nodal points re- 
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quired for the r solution and 45 for the z solution. A complete computer plot of the 
results is shown in figure 21. Two iterations were required for the z solution as well 
as for the r solution. The sum of the squares for the r solution at the final iteration 
was 0. 001. The table is given in appendix E. Again, the value of Kcp was arbitrarily 
chosen as 3.3. 

The z solution corresponding to the r solution was obtained and is also presented 
in appendix E. 


SURFACE-TENSION- DOMINATED MODEL 


As previously mentioned, the r formulation for the surface-tension-dominated flow 
can be solved independently of the z formulation (fig. 17). The initial case examined 
had a Weber number of 4, and the ratio of the radius of the jet to the radius of the disk 
was 1/2. Also, the first solution to this problem assumed only axisymmetry. The re- 
sults of the numerical solution, when a coarse mesh was used, resulted in a symmetri- 
cal r solution - symmetrical about the equipotential line emanating from point D. This 
simplified the problem in that, not only can axisymmetry be assumed, but also mirror- 
image symmetry (i. e. , symmetry about the z = 0 plane). This is significant for prob- 
lems in which the surface tension effects are to be taken into account since the surface 
tension forces are highly dependent on the curvature of the free surface. By taking ad- 
vantage of the symmetry involved, additional nodal points can be placed on the free sur- 
face without using extended computer storage. 

Referring to figure 17(a), if a vertical line were to be drawn (equipotential line) 
from point D, where r = RD, the intersection of this equipotential line with the free 
surface AG would be defined as point M. Along DM , it is known that z = 0. The vari- 
ation of r with i f/ can be computed along this line. The equipotential line DM is 
depicted in figure 22. Since DM is an equipotential line, the velocity along this line 
must be equal to V, the incoming jet velocity, with the exception of the point r = L in 
the physical plane. At point D, a velocity discontinuity will exist. 

In the physical plane, on DM , 


u = 0 


v = -if) 


on z = 0, L < r < R 


max 


(162) 


where R max is the maximum radius of the liquid flow pattern. 

An expression for the stream function along DM can now be derived, since the 
radial velocity component along DM is 0: 
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(163) 


- 1 ^ = 0 
r 9z 

This implies that \p = \J/(r) alone. To find what the function is, the definition of the 
axial velocity is used, namely, 


v =Ai£ 

r 9r 


On DM, 


y/_ 1 # 

r 9r 


Integrating this gives 


xfs =- 



+ Constant 


(164) 


(165) 


By applying the boundary condition that \]s = 0 at r = L, we can calculate the constant 
in equation (165). The following expression results: 


* =- 


yr - 2 t y l 2 

2 2 


(166) 


If this equation is nondimensionalized, 






(167) 


Dropping the starred notation, as was done in the past derivations, gives 




(168) 


Solving for r yields 
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(169) 


r = + 


^2\p + (RD) 


2 


where 



= RD 


(170) 


Equation (169) represents the boundary condition used along DM in the inverse plane 
for the r formulation. The condition z = 0 was used for the z formulation. 

In the course of finding the solution to the z formulation (fig. 17(b)), it was neces- 
sary to solve a cubic equation for the parameter T along the free surface. Physical 
as well as mathematical interpretation was required when choosing the proper root of 
the cubic since a possibility of three real roots existed. In equation (D78), appendix D, 
the only mathematically meaningful roots were those in which the absolute value of the 
parameter T was greater than or equal to 1 . However, the possibility still existed 
that all three roots would be real and in addition, satisfy the requirement that their ab- 
solute values be greater than 1. As a result, some physical insight was required when 
deciding which roots to use in the equation relating T to the fictitious point f 
(eq. (D84)). For example, it is known that, as the nozzle exit is approached along the 
free surface, z = f(r) becomes steeper (i. e. , f'(r) approaches infinity). This would 
coincide with the curvature terms dropping out of the boundary condition in the physical- 
plane formulation. An alternative approach to viewing this is that the fictitious point f 
approaches the image point z^. As a result, the value of the parameter T approaches 
1. The algorithm selected for choosing the proper root of the cubic was to select the 
value of T closest to 1 but to ensure that its absolute value was greater than or equal 
to 1. During the course of finding the solution, problems in implementing this algorithm 
occurred, particularly close to the nozzle, since the T values closest to 1 were 
slightly less than 1 and were automatically discarded by the algorithm. The resulting 
potential lines and streamlines appeared to be inaccurate when plotted. The only way 
found to circumvent this problem was to set up an additional algorithm that set T iden- 
tically equal to 1 for several free-surface nodal points in the vicinity of the nozzle (i. e. , 
those nodal points in which the computed r value was ^ 1. 009). Physically, this 
reasoning is justified since it is known that f f (r) must approach infinity there. 

The results of the numerical solution are shown in figure 23. The Weber number 
for this solution was 4, and the ratio of the radius of the liquid jet to the radius of the 
disk was 1/2. A value of 3.5 was chosen for K(p, and D was set equal to 0. 0204. In 
appendix E, the computed r and z values for the nodal points are given. 

_ P 

The sum of the squares for the r solution was 0. 77x10" , and since the r solution 
was independent of the z solution, a second iteration was unnecessary. 
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DISCUSSION OF RESULTS 


The numerical solutions were compared with the available semianalytical results of 
Schach (ref. 29) for the case of an infinite flat plate. The method employed by Schach 
is attributed to Trefftz. Schach’s method did not appear to be readily extendable to 
more- complicated geometrical flows and could not account for the effects of surface 
tension. In making the comparison between reference 29 and the numerical results, the 
fine-mesh solution presented in figure 19 was used. The results of the comparison for 
an infinite flat plate are shown in figure 24. The data points in the figure were obtained 
from figure 11 of Schach’s paper by using an expandable scale. As a result, there is 
some unknown error associated with the process of taking the results from the reference. 
In any case, the agreement looks particularly good, with the sole exception of the first 
r coordinate greater than 1. 

One final point with respect to the infinite-flat-plate solution concerns the extreme 

left coordinate in figure 24, namely, r = 4, z = 0. 125. These two values are fixed by 

continuity, both in this numerical program and in the semianalytical results of Schach. 

This result was obtained as follows: assuming constant density, the volumetric flow 

rate into the control volume at AB must balance the flow out of the control volume at 

2 

GE . In physical coordinates, the flow in is given by jtRqV and the flow out by 
27rRj e j.ZgV. Equating these and canceling lead to the fact that z^ must equal R^/2Rj e ^.. 
Nondimensionalizing with respect to Rq yields 


z* = 
g 2Rf. 


jet 


(171) 


Dropping the starred notation and observing that Rj^ = 4 at the left boundary of the 
control volume show that Zq must equal 1/8. 

As far as the finite plate is concerned, there were no available comparisons with 
past experimental or analytical work. As a result, comparisons were made with respect 
to the zero-gravity experimental data obtained in this work. The results are shown in 
dimensionless coordinates in figures 25 and 26. Figure 25 is for the case where the 
ratio of the jet radius to the disk radius is 1/2, and figure 26 for a ratio of 3/4. The 
comparisons were made with respect to the outer (or top) free surface since it was im- 
possible to view the lower free surface because of the way in which the flow occurred. 
The results were generally good for both ratios compared. Experimental data points 
were obtained from both sides of the axisymmetric sheet as it flowed around the disk. 

The average values to the left and right of the jet centerline were used to plot the solid 
lines in figures 25 and 26. The analysis corroborated the experiments in the sense that, 
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as the ratio R Q /L becomes smaller, the jet appears to leave the disk in a more nearly 
tangent manner (fig. 27). 


SUMMARY OF RESULTS 


An experimental and analytical investigation was conducted to determine the free- 
surface shapes of circular liquid jets impinging normal to sharp-edge disks in zero 
gravity. The test liquids were anhydrous ethanol and trichlorotrifluoroethane (Freon 
TF). Nozzle radii were varied from 0. 25 to 0. 75 centimeter, and disk radii of 1. 0 and 
1.5 centimeters were used. The jet velocity was varied between 12 and 365 centimeters 
per second. Under the stipulation that the nozzle be located at least 5 centimeters from 
the disk, the investigation yielded the following results: 

1. It was analytically determined that flow regions exist where viscous forces are 
not significant in the computation of free- surface shapes. It was shown that the Reynolds 
number pVRq/ ju and the ratio of jet to disk radii Rq/L uniquely define the flow re- 
gions. It was further shown that the Reynolds number specifying the transition between 
viscous and nonviscous flow decreased with an increasing ratio of jet radius to disk 
radius . 

2. Within the inviscid region, three distinct flow regimes were experimentally 

p 

found that depend uniquely on the Weber number pV Rq/<x and the ratio of jet to disk 
radii Rq/L. These flows were defined as surface tension ilow, transition flow, and 
inertial flow. The critical Weber number between regimes also decreased with an in- 
creasing ratio of jet to disk radii. 

3. A numerical solution yielding free-surface shapes and streamlines was obtained 
for the case of impingement normal to an infinite flat plate and compared favorably with 
semianalytical techniques in the literature. 

4. A numerical solution yielding free-surface shapes and streamlines was obtained 
for inertially dominated flows at ratios of jet to disk radii of 1/2 and 3/4. The compar- 
ison with experiments showed good agreement for the upper free surface. 

5. A surface-tension-dominated flow was formulated and solved numerically. The 
system Weber number was 4, and the ratio of jet to disk radii was 1/2. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, May 6, 1977, 

506-21. 


54 



APPENDIX A 


DETAILED CALCULATIONS OF FREE-SURFACE BOUNDARY CONDITIONS 


Two boundary conditions are required to solve for a free surface in a fluid dynamics 
problem. In comparison, for known boundaries, only one boundary condition is required. 
The two conditions to be satisfied are 

(1) Conservation of energy along a streamline 

(2) Zero velocity normal to the streamline 
Consider the following geometry: 


i 


I 




where i represents any point on the free surface z g = f(r) and d designates the refer- 
ence point, which is chosen as the point where the liquid jet exits from the nozzle. 


CONSERVATION OF ENERGY ALONG A STREAMLINE 

Bernoulli's equation written between points i and d becomes for inviscid incom- 
pressible, steady flow. 


— (u 2 + v 2 ) + — = ~ + — on z = f(r ) (Al) 

2 p 2 p s 

The pressure at point i, p^, is not the same as the pressure at point d, p^, because of 
the effects of surface tension. In general, 

p g " p = aJ ( A2 ) 

where p is the known gas pressure and 

o 
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J = — + — 
R 1 R 2 


where R^ and Rg are radii of curvature and 


_1 

R 


1 / 

1 (! + f 2 ) 


R 2 / *2 
z r(l + f J ) 


Combining equations (A4) and (A5) with (A2) gives 


(A3) 


(A 4) 


(A5) 


P g - P = a 


(l + H) r(l + f 2 ) 


(A6) 


which can be combined to yield 


ad/ rf' 
Per “ P = - — 


^ r dr 


vVT 


(A7) 


+ f 


If at point d we apply z = H, r = Rq, and df/dr = f' = 00 , then for large f', ^1 + f ^ = 

^2 

= f*. Therefore, substitution into equation (A7) yields 


p e- p d= — 


R, 


(A8) 


Similarly, at point i, 


ad / rf' 

p g"Pi = - — f 


r dr 




+ f 


(A9) 


Now equations (A8) and (A9) can be substituted into equation (Al) to obtain 
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a 


(A 10) 


1 ( U 2 + V 2).^_dj 

2 pr dr 


rf’ 


<Vi - f ' 2 y 


= 1^2 . _2_ 

2 pR 


on z g = f(r ) 


0 


ZERO VELOCITY NORMAL TO FREE SURFACE - 

A portion of the free surface z = f(r) is shown in sketch (b) 

s 



■ f(rl 


where n is the unit outward normal to the surface at some point and u and v are the 
velocity components such that V = ui + vj . Since the velocity normal to the surface is 
zero, 

V • n = 0 (All) 

With the surface given by z g = f(r), the unit normal is given as 

n = - -*!=-!_ (A 12) 

fi T7 2 


Taking the gradient at the surface and dividing that result by the absolute value of the 
gradient (i. e. , h = V<p | Vcp | , where i is a unit vector in the radial direction and j is 
a unit vector in the axial direction, yield 


V • n 


uf* 

VT77 2 


+ 



(A13) 


which simplifies to 


-uf' + v = 0 on z g = f(r) 


(A14) 
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APPENDIX B 


ZERO-GRAVITY DROP TOWER TEST FACILITY 

The experimental data for this study were obtained in the Lewis Research Center's 
2.2-second zero-gravity facility. A schematic diagram of this facility is shown in fig- 
ure 28. The facility consists of a building 6. 4 meters square by 30. 5 meters tall. Con- 
tained within the building is a drop area 27 meters high with a cross section 1.5 meters 
by 2. 75 meters. 

The service building has a shop and service area, a calibration room, and a 
controlled-environment room. Those components of the experiment that required spe- 
cial handling were prepared in the controlled-environment room. This air-conditioned 
and filtered room contains an ultrasonic cleaning system and the laboratory equipment 
necessary for handling test liquids. 


MODE OF OPERATION 

A 2. 2-second period of weightlessness is obtained by allowing the experiment pack- 
age to free fall from the top of the drop area. To minimize drag on the experiment 
package, it is enclosed in a drag shield designed with a high ratio of weight to frontal 
area and a low drag coefficient. The relative motion of the experiment package with 
respect to the drag shield during a test is shown in figure 29. Throughout the test, the 
experiment package and drag shield fall freely and independently of each other; that is, 
no guide wires or electrical lines, etc. , are connected to either. Therefore, the only 
force acting on the freely falling experiment package is the air drag associated with the 
relative motion of the package within the enclosure of the drag shield. This air drag 
results in an equivalent gravitational acceleration acting on the experiment, which is 
estimated to be below 10”^ g's. 


RELEASE SYSTEM 

The experiment package, installed within the drag shield, is suspended at the top of 
the drop area by means of a highly stressed music wire attached to the release system. 
This release system consists of a double-acting air cylinder with a hard-steel knife edge 
attached to the piston. Pressurization of the air cylinder drives the knife edge against 
the wire, which is backed by an anvil. The resulting notch causes the wire to fail, 
smoothly releasing the experiment. No measurable disturbances are imparted to the 
package by this release procedure. 
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RECOVERY SYSTEM 


After the experiment package and drag shield have traversed the total length of the 
drop area, they are recovered by a deceleration in a 2.2-meter-deep container filled 
with sand. The deceleration rate {averaging 15 g's) is controlled by selectively varying 
the tips of the deceleration spikes mounted on the bottom of the drag shield (fig. 28). At 
the time of impact of the drag shield in the deceleration container, the experiment pack- 
age has traversed the vertical distance within the drag shield (compare figs. 29(a) and 

(c)). 
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APPENDIX C 


DERIVATION OF INVERSE- PLANE EQU ATIONS 
FORMULATION 


Assume that there exist functions ip = ip(r, z) and <p = <p(r,z). Then, there are 
inverse functions r = r (ip, ip) and z = z (<p, ip) such that 


9r 1 dip 

d(p J dz 

9z _ _ 1 9 (£_ 

dip J 9r 

9z _ 1 dtp 
dip J 9r 


9r _ _1 9 <p_ 

dip J 9z 


and where the Jacobian 


J is defined as 


dip 

dip 

dr 

dz 

dip 

dip 

dr 

dz 



(Cl) 

(C2) 

(C3) 

(C4) 


(C5) 


From these equations follow two important relations. Substituting the values of dcp/dr 
and dip/dz from equations (Cl) and (C2) into equation (112) yields 


9z _ _ 1 _9r_ 

dip r 9 <p 


(C6) 


Similarly, substituting the values for dcp/ 9z and dip /dr from equations (C3)and (C4) 
into equation (113) yields 


9r _ 1 9z 
dip r d<p 


(C7) 
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For subsequent relations, it is necessary to express the Jacobian in terms of the 
inverse functions r {ip , cp) and z(\p, cp). Substituting equations (Cl) to (C4) into equa- 
tion (C5) results in 


J = j 2 \(J*\ (Jl\ _ /Jr \ / 9z \1 
J_\a ify/Xdcp/ \dpl\dtp/\ 

Now, using equations (C6) and (C7) in equation (C8) yields 


J = - 



(C8) 


(C9) 


GOVERNING EQUATIONS FOR INVERSE FUNCTIONS 


Differentiating equation (C6) with respect to p yields 


a 2 z 


d\p 2 r d ^ d( P r 2 8 & 8< P 


1 3 2 r + 1 3r 3r 


(CIO) 


and combining this with the derivative of equation (C7) with respect to cp, which is 


3 2 r 

d cp d\fy 


= - 1. (*_)(*)+ k 

2 \3c o) \3 wj r 


a 2„ 
3 z 


r _ Kdtp/ \dcpt r d(p 2 


(Cll) 


gives 


3 2 z 


1 3r 3z 


3 2 z + _1 


3 p 2 r 3 d<p d(p r 2 3 <p 2 


_ Jr_ j£_ 

2 dtp dip 


(C12) 


By using equations (C6) and (C7), the following important equation is obtained: 


r 2 d 2z + 9 2 z 2 3z 3z 


3i// 2 3<p 2 


d\f/ dcp 


(C13) 
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On the other hand, differentiating equation (C6) with respect to cp and equation (C7) 
with respect to and combining the result lead to the equation 


0 2 r = _ j?r_ + J_ /5r\ 2 . _1_ _0r _0z_ , { 

d\[/ 2 r 2 d(p Z r 3 ' 9<59 ' r 2 9 ^ 8<p 

By substituting for dz/dcp from equation (C7), the following equation for r(cp , \p ) is 
obtained: 


S 2 r + 1^ _ J_ / 9r \ 2 + 1 / Sr \ 2 _ Q 

d * 2 r 2 d<p 2 r 3 ^/ rW/ 


(C15) 
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[ 


APPENDIX D 

DETAILED DERIVATIONS OF FINITE DIFFERENCE 
OPERATORS ALONG BOUNDARIES 

INFINITE FLAT PLATE 

The finite difference representations for the derivatives are those shown in the text 
(eqs. (142) to (149)). Refer to figure 13. 

— = 0 on Line EC 

3 \p 

The line EC is shown in sketch c. Applying the r difference equation (eq. (150)) 



2 

3 

0 1 

<b 

, & 


(0 


at point 0, where f is a fictitious point outside the boundary, gives 
.3 


r 0 " ( r 2 + r f) + r 0 


— " ~ <r 2 - r / 


a 


2 8 


-5- (r x + r 3 ) + -J— (rj - r 3 ) 2 = 0 (Dl) 


2a" 


8a" 


Along EC, 3r/3 ip = 0. This implies that 


r 2 ~ r f 
2 A <// 


= 0 


(D2) 


Therefore, 


r 2" r f 


(D3) 


63 


I 


Equation (Dl) then becomes 


r 


4 

0 _ 


rXr 


0 2 



l^( r l + r 3 ) + ^< r l- r 3 ) 2 = 0 


2 a‘ 


8a 


(D4) 


— = 0 on Line GE 

dxf/ 


Applying dr/dxp on GE (sketch d) yields *2 =r 4 anc * < ^ oes n °t involve the unknown 


I© 

2 

f 1 

0 " 

4 

te 

(d) 


fictitious point f. However, dr/d\p =0 implies that r*r (\p). Therefore, r must 
only be a function of cp, r =r (<p). But, along GE, cp is constant (cp =0). Hence, r 
is constant along GE. 



on Line AG 


For line AG (sketch e), the r difference equation can be written 


© ? 

f 

0 1 © 


,4 


(e) 
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1 


_4 0 / \ 2 

r 0 ‘ Y f 4 0 


-L-±(r f - r„) 2 
a 2 8 f 4 


< r l + r 3> + < r l - r 3> 2 
2 a z 8a * 


Along AG, 


3r\ 2 + r 2 /Jr f „ , 



0<p/ \9^ 



This implies that 


i - r 3\ 2 . _2/V r 4\ 2 


2 A<p 


+ r 


0 


2 A^/ 


= 1 


Rearranging and letting Ai // = a (recalling that a = Acp/Axp ) gives 


( r f “ r 4 )2 = — 

r 2 
r 0 


~ < r l - r 3) 2 

a 


Upon which we can obtain 


r f - r 4 = 


- lA 2 - ~ 

r 0 I a 2 


( r l - r 3> 


or, rewriting, 


r f + r 4 = 


— V 4 * 2 - i 

r 0 y a 2 


(r x - r 3 ) + 2 r< 


Inserting equations (D8) and (D9) into (D5) yields the desired relation 


4 3 2 

r 0 ' r 0 r 4 + r 


0 — - I / 4 a 2 - J- 

a 2 2 1 a 2 


(r x - r 3 ) 2 2 . (r, + r,) 


1^ 


2a‘ 


+ i 

8 


— (r t - r 3 ) 2 - 4a 2 


a 


0 (D5) 


(D6) 


(D7) 


(D8) 


(D9) 


(DIO) 
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Constant z on Line AB 


For z to be constant on AB , let z = K <p and dz/dcp =0 on BC (sketch f). Apply- 


© 





(f) 


ing the z difference equation (eq. (151)) at point 0, where f is a fictitious point outside 
the boundary, gives 

Zf)( Z l - Zg) 

(Dll) 

Along BC, 3z/9\^ = 0. This implies that z 2 - z f = 0 or 


z n = 



+ z f ) + 


2(r lot 2 + 


(z 


1) 


1 + z 3 


) + 


( r »“ + i) 


( z s 


Zr — Z c 


(D12) 


Therefore, equation (Dll) becomes 


z 0 = 


r 0 z 2 


r 2 + J_ 2(r 2 oa 2 + l) 


(Zj + Zg) 


(D13) 


= 0 on Line GE 

d cp 

For the line GE (sketch g), the z difference equation can be written 
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© 


2 

f 0 1 


4 


;© 




& 

9 -J\ ( z 2 + z 4> + — < Z 1 + z f> + ~jY A (Z 2 " Z 4 )(Z 1 " Z f } 


2(rja* + 1) 


( r »" + i) 


(D14) 


Along GE, dz/dcp = 0. This implies that Zj - z f = 0 or 


z f sz l 


(D15) 


Therefore, equation (D14) becomes 



2 2 , 
r 0 a + 1 


(D16) 



on Line AG 


For the line AG (sketch h), the z difference equation can be written 


©3 


,f 

0 


(h) 
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z o = 


fo'-h 


0 ( Zf +Z 4 ) + -1 (Zj + ZgU— — - ^( Z f Z 


2(r?a 2 + 1) 


( r »* + i) 


Along AG, 


Hence , 


Y— f +(— f = 

\9i/// \dcpl 


r 2/Vf4\ 2 7 Z| - Z 3\ 2 = i 

U \2 A i// / \2 A (p 


which can be expressed as 


(z f - z 4 ) 2 =i- 


0 L 


4a 2 - — (Zj - z 3 ) 2 


a 


This yields 


z f - z 4 =. 


-JR 71 

r o yL *2 


( Z 1 - z 3^ 


z f + z 4 - 


i 


4a 2 - — (z j - z 3 ) 2 

O' 2 


+ Z, 


Inserting these last two equations in (D17) eliminates the fictitious point z 


z 0 = 


0 


V'i) 


A. J4, 2 - -L 

r 0 V a 2 


(Zf - z 3 ) + 2z^ 


2(r 2 a 2 + 1) 


(Zi+Zg) 


+ - Z - r Z 3 . L 2 --! 

4r 0 ( r 0 « + V « 2 


( z 1 


4 )( z l - z 3 ) 

(D17) 


(D18) 


(D19) 


(D20) 

(D21) 

f and yields 


z 3 ) 2 (D22) 
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A special boundary condition is required for point G since it is a part of two 
separate boundaries (sketch i). Applying the z difference equation yields 


f 2 ©T z i 


© 


(i) 


© 



+ z 4 ) + 


.2, 

0 C 


2(r ^aT + 


1) 


( Z 1 + f 2 


) + 



z 4^ z l " f 2> 


(D23) 


Applying the boundary condition along GE , namely, 3z/3 <p = 0, gives fg = Zj. There- 
fore, 


z 0 = 


2 f» + -l 


0 (fj + z 4 ) + 1 


2 2 i 

r^a + 1 


The boundary condition along AG (eq. (D20)) gives 



but f z = Zj. Hence, 


f 


1 


z 4 + 


2a 

r 0 


(D24) 


(D25) 


(D26) 
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Equation (D23) becomes 


r 0 z 4 + 


z 0 ~ 


2 1 2 2 , 
r o + — r 0 O! +1 

or 


(D27) 


such that Zq = Zq(zj, z 4 ) at point G. 

FINITE PLATE (EXCLUDING SURFACE TENSION) 

The differences in the boundary conditions between the infinite and finite plates 
involve GE and the addition of the free surface ED. In addition, G becomes a special 
point in the r formulation, and both G and E become special points in the z formu- 
lation. 



= 0 on Line GE 


Along the line GE (sketch j), the r difference equation is 


t© 

2 

l 0 _ 1 

4 

m 


r 0 - J (r 2 + r 4> + r 0 


— - - (r 2 - r 4 ) 2 
c 2 8 . 


2 or 


(r 1 + f) + J_( ri - f) 2 =0 (D28) 
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In finite difference form, the boundary condition, with 


dz _ z 2 ~ z 4 
dr r 2 - r 4 


is 


r 1 - f r n 

- 1 + - 2 -( z 2 - z 4 )=0 

A (p A\f/ z % 


(D29) 


In equation (D29), central finite difference was not used. Solving for f gives 

f =r l + ar o( z 2 ' z 4> 

and inserting this into the r difference equation yields the desired result: 


r 0 " J (r 2 + r 4> + r 0 




r °- [2r t + o-r 0 (z 2 - z 4 )] + ^ (z 2 - z/ = 0 

(D31) 


J 2a 


2 2 

( + r^ ( — \ = 1 on the Free Surface ED 
\a w) \dxl// 


On the free surface ED (sketch k), the r difference equation can be written 


© 


i ® 


(k) 


r 0'y (r 2* r f )+r 0 


Xi (r *' r * )2 


2a i 


(r 1 + r 3 ) + — — (r 1 - r a ) = 0 


8a 


1 A 3' 


(D32) 


71 



Using 


on ED gives 



which yields (with a = A \p) 


( r 2 - r f ) 2 



(D33) 


(D34) 


and 



Substituting equations (C34) and (C35) into (C32) gives 


(D35) 


4 3 2 

r 0 " r 0 r 2 + r 0 




( r l + r 3 ) 






(D36) 


Now, points E and G will be special points since equation (D31) cannot be directly 
applied there. One of these positions can be specified as known, r^ = Constant. Let us 
examine the special point G (sketch l). Applying equation (DIO) at point G yields 
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I 


m 

© 


id 




3 2 

- r Q r 4 + r 


0 





( r l + f 2> 



A (r 1 - f 2 )2 - 4a 2 


= 0 


Now, applying 


/A\ + r (AL\/dz\ = o 

\d<pj \3^//\dr/ 
along GE without involving f ^ gives 

Note again that central finite difference was not used to derive this expression, 
for fg gives 


f 2 =r 1+ 2ar 0 (z 0 - z 4 ) 
Therefore, equation (D37) becomes 


r 


4 

0 


r 0 r 4 


+ r, 




*aY 


r o r , 

" L r i + ar o( z o 

OL 



+ 


r 0 (z 0 




(D37) 

(D38) 

Solving 

(D39) 

0 (D40) 
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r /J2L\ _ /9£\/dz\ = o Qn Line GE 

Xdxjy/ \d(p/\<3r/ 


Along the line GE (sketch m), applying the z difference equation yields 


i© 



■VSl/ 

2 

f 

0 1 


4 


© 


(m) 


z 0 = 


0 (z 2 + z 4 ) + Zl + * . < z 2 - z 4>< z l - f > 


2 fo + -l 


2(rQO! 2 + 1) 4 (v 2 Q a + 


Now, along GE, 


\3i/// \d(p)\dr) 


which in difference form can be written 


< z 2 ~ z 4' 


' Z 1 - { \/ z 2 ' z 4' 


2 Ai^ / \2 A.<p /\r f 


= 0 


Solving for f gives 


f = z i ~ r 0 a < r 2 - r 4 } 


(D41) 


(D42) 


(D43) 
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Hence, equation (D41) becomes 


z 0 


2 To + ^ 


( z 2 + 24 ) + 


2z l " r o<*( r 2 - r 4 ) | ( z 2 - z 4 } r 0“ (r 2 ' r 4> 


2(r 2 « 2 + 1) 


a 


('!*•;) 


(D44) 


2 2 

r 2 ^5-^ + (-^\ = 1 on Free Surface ED 



d\p/ \3 tp 



On the free surface ED (sketch n), the z difference equation can be written 


© 3 

2 

0 1 

© 


>f 


(n) 


Z 0 = 




( z 2 + z f) i 7 ( z i + z ^) ~ r ( z 9 ” z f)( z i “ z ^) 


2(r 2 <* 2 + 1) 


J 1 T 3 


( r °“ + t) 


J 2 f A 1 ""S' 


(D45) 


Applying 


on ED gives 


r 2 /Sz\ 2 + /3z f = J 


3^/ W 



r 


2 

0 




= 1 


(D46) 
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Substituting equations (D47) and (D48) into (D45) yields the desired expression 




Applying equation (D22 ) at point G (sketch o) yields 
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4a 2 - — (z 1 - f 2 ) 2 + 2 z 4 + 1 — 


2(r 2 n a 2 +-1) 


Z 1 I f 2 / 4a 2 . _1 

: 0 po«a)V - 


— — ( z i - f 2 ) (D50) 

2 l 2 


To find f 2 , we apply 


r _ /_az_\/dz\ 0 


d\}// \d(p/\drS 


Without involving f^, 




(D51) 


Solving for f 2 , 


f 2 =Z 1 " 2r 0 a(r 0- r 4 } 


(D52) 


Substitution into equation (D50) yields 


r 0 1 J 2 2, ,2 Z 1 - r o“< r O - r 4 

- ya -r 0 (r 0 -r 4 ) + z 4 + — 

r 0 + "T r 0 J r 0“ 2+1 

Q! 


' 0 " r 4 */„2 „ 2 ,„ „ ^2 

- ya - r 0 (r Q -r 4 ) 


(D53) 


Let us examine the special point at E (sketch p). Applying equation (D49) at point E 


77 


results in 



© 


2 

[z 0 

© 1 ® 


f l 

(p) 



i /^7X 
r 0 V c 2 


(Zi - f 9 ) 


2(r 2 a 2 + 


1) 


(Zj + f 2 ) 



In order to find the fictitious point fg, we apply 


r 





= 0 


on GE, without involving f^: 

r p ( z 2 ~ z q ) _ / Z 1 ~ f 2\ / z 2 ~ z 0 \ = 0 

\2A?/\r 2 - r J 


(D55) 


Solving for fg, 


f 2 =z l“ 2r 0 a < r 2 “ r 0 ) 


(D56) 
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z 0 = 


-f V a2 - r 0< r 2 - 


tj' 


Z 1 - r Q tt ( r 2 ~ r 0> 

' 2 c 2 
O c 


aT 


r 0 + 




a 2 - r 2 (r 2 - r Q ) 2 


(D57) 


a 


FINITE PLATE (INCLUDING SURFACE TENSION) 



1 


1 - 


2,21c! 
We We r dr 



on Line AG 


Applying the r difference equation (eq. (115)) at point o on the line AG (sketch q), 




f 

0 


I 1 

<q) 




where f is a fictitious point outside the boundary, gives 



— - ( r 
2 


f + r 4 


) + r 


0 



-r(r# 




~^±(r 1 + r 3 ) + -1- (rj - r 3 ) 2 = 0 (D58) 
2a 2 8a 2 


In difference form, the free-surface boundary condition becomes 
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1 


(D59) 



1 - 


2 | 2 1 d 

We We Tq drQ 



i-r 3 )/2A <p\ 2 

+ r 0 


,( r f - r 4 )/2 Axp j 


Simplifying by using a = Acp/A\p and defining 

' < r l - r 3> 

Q = ±2— 

(r f ~ r 4> 

yields 


“ ( r l " r 3 ) 2 + r o( r f " r J 2 =' 

or 


4a" 


1 - 


2 1 d 


We 


We r Q dr Q 




(D60) 


Now, examine 


r O 



As an approximation to this derivative Q is assumed as a constant. In actuality, 
Q = f(rj) and r^ = f(rQ). Expanding 


1 d 
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yields 


1 d_ 

r n dr. 


2Q + r 2 


° ° \VQ + r o/ V Q + i*o (Q + r J) 

Therefore, the finite difference representation for the free surface becomes 


1 /„ v2 2/ \2 _ 

— - (r x - r g ) + r Q (r f - r 4 ) - 

Cl 


4a 2 

1--L + J- 

2Q + rg 

We We 

J)/q + Tq (Q + r|) 


The term must be eliminated between equations (D58) and (D62). If we define 


x = r f - r 4 


equation (D62) can be written 


1 , \2 2 2 

— (r, - r 3 ) + = 

a 2 


4a 


1- 2 + 2 

4a 2 

We We 

|/q + r 2 (Q + r 2 ) 


where 


Q = 


( r l - r 3 ) 


(?) 


Inserting (D63) into (D58) yields 


2 2 3 

r 0 x r 0 


.4 . _3. 


a + “ x - r o + r 0 r 4 - + (r l + r 3 ) - < r i - r 3> 2 = 0 

s 2 a 2 2a 2 8a 2 


(D61) 


(D62) 


x as 
(D63) 


(D64) 


(D65) 


(D66) 
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This is of the form 


A’x 2 + B'x + C* 


= 0 


Hence, 


X 


B’ + 


2 A 1 


4A*C' 


where 


A' = — - 
8 


r 3 

B*=!® 


c* = ” r Q + r 3 ^ 


— — + — — ( r i + r„) — — (r , - r»f 


a 


2 a 


8 a 


(D67) 


(D68) 


(D69) 


(D70) 


The sign in front of the square root in equation (D67) was chosen as positive since x 
must be greater than or equal to zero. In addition, employing equation (141) and the 
results of appendix C yield 



1 


1 2 ( 2 1 d 

r 

We We r dr 

J r 2 dz/d* 2 + j 
_ f dz/dq> 


on AG . 

Applying the z difference equation (eq. (114)) at point 0 (sketch q), where f is a 
fictitious point outside the boundary, yields 
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z 0 - 


2 ( r 0 + 


^ JT (f + z 4 ) + 


(f - Z 4 )(Zj - Zg) 

(Zj + Zg) + 

2(r 2 a 2 +l) 4 1 


a 




In finite difference form, the free-surface boundary condition can be written 


rj)(f - z 4 ) 2 + (zj - z 3 ) 2 - 


4a 


a 


1 - 2 | 2 1 d 


We We r n dr 


odr °\f 2 o- 


IQ* + 


where 


Q* = a- 


2 ( f " z 4^ 

( Z 1 " V 


As an approximation Q* is assumed to be a constant. Expanding 


1 d 
r 


0 


° dr ° W 


rgQ* + iy 


gives 


1 d 


^ dr ° V r 0 Q * + X j r 0 (r 0 Q * + !) V r 0 Q * + 1 


Therefore, the finite difference representation along the free surface becomes 


r 0 (f " z 4 )2 + \ (z l " z 3 )2 " 2 2 1 

o- z 1 _ z 1 


4a 


1 


We We r 


0 (r 2 Q* + 1) 


V r 0 Q* + 


(D71) 


(D72) 


(D73) 


(D74) 


(D75) 
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The point f must now be eliminated between equations (D76) and (D75). This is done 
as follows: Equation (D75) is solved for f = fct(zj, Zg, z^). The results are then in- 
serted into equation (D71). In this way, equation (D71) remains explicit in Zq. Ac- 
tually, it will be more convenient to solve for the variable (f - z^) instead of f since 
equation (D71) can be expressed as 


z 0 = 


r 0 z 4 


z l + z 3 


r 2 + J_ 2 (r 2 0 « 2 + l) 

or 


(f - z 4 ) 


2ari+ (z 1 - Zn ) I 


4 a 




(D76) 


Equation (D75) is solved for (f - z^) 


2 2 , t \2 
a r 0 (f - z 4> , , 

(z, - z,) 2 


„ 2 2 
4a a 

(zi - z ,) 5 


1 - + 

We We r rt 


(D77) 


.2.2 


r o tt ( f - z 4>* 


(z, - z„r 


/ 2 2/j, x2 

r n a (f - z.) 


■+ 1 


(Zi - Zn) 


A new variable, T , is introduced: 


T 


2 


2 2/. N 2 

<* r p( f - z 4> 1 

( Z 1 " z 3^ 2 


Also, let 



(D78) 


(D79) 


(D80) 


(D81) 
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Making these substitutions into equation (D77) results in a cubic equation for T: 


I 


T 


3 



(D82) 


Rewriting this gives 


T 3 + aT + b = 0 (D83) 

Once T is found from equation (D83), (f - z^) is calculated from equation (D78) as 
follows : 


f - 



(D84) 


The choice of sign in equation (D84) depends on whether z is less than, equal to, or 
greater than zero, since the fictitious point f by definition must lie outside the bound- 
ary. 
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APPENDIX E 


COMPUTER SOLUTIONS 
INFINITE FLAT PLATE 
Coarse- Mesh Solution for r and z 


t 

G 


4 

3.659 

^ 2. 616 

°0. 125 

Zj-0. 137 

z 2 =0. 191 

4 

3.655 

2.616 

' z g =0. 100 

z 9 =0. 109 

z 10 =0. 152 

4 

3.652 

2. 616 

"z 16 =0.075 

Zjy = 0. 082 

P — ' 

OO 

II 

o 

M 

£ 

4 

3.650 

2.616 

U z 2 4=0.05 

' z 2 5=0. 055 ’ 

z 2 g=0. 076 

4 

3.649 

£616 

°Z32=0.025 ' 

^33=0.027 ' 

’234=0. 038 

4 

3.648 

2.616 

0 

E 

’o 

0 


Jl. 258 
T z 3 = 0 . 226 

^1. 525 
z 4 =0. 367 

^1- 054 
z 5 =0. 878 

l. 232 

|z u =0. 181 

1.492 
z 12 =0. 306 

0.965 
z 13 =0. 814 

^ 2. 212 
Jz 19 =0. 136 

1.460 
Zgj-a 239 

0.848 
z 21 =0. 733 

I 2 . 198 
’^27=0.091 

;.433 , 

z 2 g=0. 166 

0.702 

2^=0.625 

I 2 . 189 

Tz 35 =0.045 ' 

1.412 
^ 3 ^= 0 . 085 

0.494 

237 = 0.45 

12. 187 

.1.404 

0 

^ 0 “ 

0 

0 


C 

(r) 


A 


^ 1. 024 

1. 007 

.1-0 

V 1 - 45 

z 7 =2.35 

Kp=3. 315 

J.918 , 

0.904 

0. 895 

z 14 =1.44 

*15* » 

K<p=3. 315 

0.8006 

0. 986 

0.775 

z 22 =L43 

z 2 3=2. 34 

W=3. 315 

0.655 , 

z 30 =L41 

0.643 
z 31=2. 34 

0.632 
K<p=3. 315 

0.456 

z 38 =L36 

0.4474 

z 39 =£ 32 4 

0.447 
K<p=\ 315 

0 

z 40 = l* 10 ~ 

L o 

z 4 i=2.21 

a. 00 

K^=3. 315 


B 
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£ |( 3. 858 n 3. 690 ) ( 3. 128 ,, 2.637 )( 2.451 ^ 2. 253 )( 1.879 ^ 1.526 ^ 1. 268 Oil ^ 0.991 0.971 ^ 0.963 ,,0.955 0. 951 0. 949 

z 16 -0.H2 z 17 -0.116 z 18 -0.I22 z 19 -0. 143 z^-0. 169 z 21 -0. 182 z 22 =0.200 ' z 23 -0.Z41 " z 24 -0. 309 "z 25 -0. 454 % 26 =0. 830 29 *z 28 -1.79 "2^30 *z^l79 "z 31 -3. 26 * K^-3.315 



I 


z 140-°' 552 z M 2'L« z m -2.20 2,^176 Kp-3.315 





Coarse-Mesh Solution (Rq/L = 3/4) 


ill ■ o 


Zj— 1.40 

Zj.-0.887 

z 3 -0, 341 

z 4 - 0 . 2 % 

z 5 *0.597 

z 6 =L03 

Zj=l. 88 

z 8 =Z93 

Kp*3. 3 

y 1=4. 05 

y 2 -3.39 

y 3 “^ 67 

y 4 *L55 

y 5 -1.05 

V 1 - 00 5 

y ? =0.995 ' 

y^0994 ' 

r*l 

z 9 -L417 

z 10 =-a 906 

: 

z u -0.373 

; 

Z 12 -0. 238 

z 13 -0. 505 

; 

z 14 -0.948 

z 15 *l. 83 

Zi 6 -Z 86 

K<p=3. 3 

y9=4.03 

Yio * 3 - 37 

y 11 * 2.66 

y 12 *L 50 ! 

y 13 '(l 956 

y 14 =0.897 ! 

y 15=0.889 ' 

y 16*0* 888 

r=0. 895 

* 

i— J 

1 

r-H 

— 

z 18 -0.938 

z 19 — 0.404 

z 20"0- 179 

Z2i*0, 404 

z 22 =0. 842 

► 4 

z z 3 =1.77 i 

z 24 - 2 . 77 

K^=3. 3 

y 17 »4.oi 

34 

y 19=2.64 

y 20 = !• 46 

^ 21 * 0 . 845 

Y22 = 0. 776 

y 23 =0. 768 

y 2 4=0. 767 

r=0. 775 

z 25 - 1.451 

2 26"^* 965 

z^ -0.437 

z 28 =0. 119 

Z29=0, 290 

z 30 =0.702 

Z 3 J- 1. 68 

> . 

z 32 *2. 67 

Kp-3. 3 

| 

y 2 5'4 . 00 

V 2 6=3- 32 

i 

727 “2.62 

!y 28 =1.42 

'y 29=0. 710 

\ 

y 30=0. 631 

1 

Jy 31=0. 623 

y 32 * 0 . 622 

1 

r=0. 632 

Z 33 *- 1.468 

Z34-0.991 

t i 

Z 35 -O .470 

^-0.059 

^=0.159 

z 3g =0.491 

i - 


z 40 =2. 54 

K<p=3.3 

Y33 = 3. 98 

V34=3.30 

y35*2. 61 

Y36" 1 - 37 

y 37 = 0 . 533 

*38=0.436 

y 3 9=0. 4289 

y 40 =0.428 

i r=0. 447 

z 41 — 1.485 

z 42 =-1.018 

z 4 3°-0. 503 

z =0 

> 4 

z =0 

z =0 

t 4 

z 44 =l. 10 

z 4 5 = 2 . 20 

K^=3. 3 

y 4 l=3.96 

V42 =3 * 28 

y 43 =2. 59 

r=l. 333 

r=0. 238 

r -0 

y 4 4 = 0. 0 

y 45 = 0.0 

V =0 


(u) 
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y r 1. 0003 yrl-007 y 3 -1.0015 y 4 -1.0058 yg-1.018 y 6 -l. 032 y 7 -1.432 y g =2.236 y 9 =1.432 y r 1.032 y u =l. 018 y ir l. 0058 y ir 1.0015 y 14 =1.0007 y 15 -1.0003 



y I51 -0. 2471 ^y I52 =0. 4460 ^y 153 -l. 348 2 . 000 y 154 =l. 348 y 155 =0. 4460 y^O. 2471 0 . 0 







APPENDIX F 


SYMBOLS 


A 


a 


B 

b 

C 

C 0’ C l> C 2 
c 


D 

d 

d C 

E 

F 


ir 


G 


g 

H 

i 

M 

J 

Kcp 

L 

N 

n 

O 


numerical constant 

stream function increment, numerical constant 
numerical constant 
numerical constant 

numerical constant; integration constant 
numerical constants 
numerical constant 

1/a 2 

designated reference point 
nozzle diameter, cm 
special point 
function 

unknown scale factor, m 
functions 

fictitious points; functions 

special point 

gas 

distance between plate and nozzle, m 
any point on free surface 

unit vectors in radial and axial directions, respectively 

Jacobian; local curvature 

numerical constant 

disk radius, m 

index 

unit normal on free surface 
stagnation point 
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p 

00 

Q 

Q* 

R 

max 

R o 

R 1» R 2 

RD 

Re 



u 


•r 

v 

v o 

V 

We 

z 

z o 

a 

ft 

6/6r 


pressure, N/m 2 

o 

unknown scale factor, N/m 

< r l - r 3> 2 
( p f - r 4 ) 2 

« 2 (f - z 4 ) 2 

< Z 1 " z 3^ 2 

maximum radius for surface tension model, 
radius of nozzle, m 
radii of curvature, m 

L/R 0 

Reynolds number, pVRq/ju 
radial coordinate, m 
radial dimension, m 
unknown scale factor, m 
nozzle radii, m 

a 2 rQ(f - z 4 ) 2 /(z 1 - z 3 ) 2 + 1 
unknown scale factor, m/sec 
radial velocity component, m/sec 
jet velocity, m/sec 
vector velocity, m/sec 
unknown scale factor, m/sec 
axial velocity component, m/sec 

o 

Weber number, pV Rq/ct 
axial coordinate, m 
unknown scale factor, m 

A.cp/A.xjy 

3 2 

specific surface tension, a/p, m /sec 
finite difference analogy of 9/3r 



M 

P 

cr 

<P 

X 

* 


9 

liquid viscosity, g/m -sec 

O 

liquid density, kg/m 
liquid surface tension, N/ m 
velocity potential, m /sec 
r f - r 4 , m 

O 

Stoke ’s stream function, m /sec 


Subscripts: 

cr 

d 


E,G 

g 

i,k,n 

jet 

max 

s 

0 

0 , 1 , 2 , 3 , 4 


critical 

datum 

special points 
gas 

free indices 
jet 

maximum 

surface 

stagnation point, subscript for z 
nodal point locations 


Superscripts: 


differentiation with respect to r 
dimensionless quantity 
dimensionless quantity 
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TABLE I. - LIQUID PROPERTIES AT 20° C 


Liquid 

Surface 

tension, 

N/cm 

Density, 

g/cm 3 

Absolute 
viscosity, 
g/cm- sec 

Anhydrous ethanol 
Trichlorotrifluoroethane 

22.3X10" 5 

18.6 

0.789 

1.579 

1.200x10' 2 
.700 
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TABLE II. - SUMMARY OF PARAMETERS 


Test 

liquid 

Nozzle 

radius, 

R o > 

cm 

Disk 

radius, 

L, 

cm 

Ratio, 

V L 

Velocity 
of jet, 
cm/sec 

Reynolds 

number, 

PVRq/m 

Weber 
number, 
PV 2 R 0 /ct 

Flow 

category 

Freon TF 

0. 

25 

1 

0 

0. 

25 

34.3 

1934 

24.8 

S 

Freon TF 







38. 0 

2133 

30.5 

S 

Freon TF 







43.7 

2464 

40.4 


r 

Ethanol a 







77.0 

1278 

52.2 



Freon TF 







50.5 

2848 

54.0 



Freon TF 







58.8 

3316 

73.4 



Ethanol 







101.0 

1662 

90. 1 



Freon TF 







79.6 

4460 

113.2 



Ethanol 







131.3 

2160 

151.9 



Freon TF 







88.6 

5000 

161.5 



Ethanol 







149.8 

2465 

197.8 











237.7 

3911 

498 











282.2 

4642 

703 





i 






365.6 

6014 

1180 



Ethanol a 

.50 

1 

5 

.33 

26.4 

868 

12.2 

S 

Freon TF 







19.5 

1815 

16. 1 

S 

Ethanol 







37. 7 

1240 

25. 1 

S 

Freon TF 







27.4 

3090 

31.8 

T 

Ethanol 







46.0 

1513 

37.2 



Ethanol 







53.9 

1773 

51.2 



Freon TF 







37.6 

4240 

60.0 



Ethanol 







59.6 

1900 

62.0 



Freon TF 

1 






42.75 

4830 

77.6 

I 


Ethanol 

.75 



.50 

15. 4 

760 

6.2 

£ 




.75 





19. 1 

942 

9.6 





.75 





19.9 

986 

10. 5 





.50 

1.0 



26.4 

868 ' 

12.2 



Freon TF 

.50 

1. 0 



17.7 

1996 

13.2 



Ethanol 

.75 

1:5 



28.2 

1891 

21.0 



Ethanol 

.50 

1. 0 



34.9 

1148 

21. 3 

T 

Freon TF 

.50 

1. 0 



23.3 

2628 

23.0 



Ethanol 

. 75 

1 

.5 



34.4 

1697 

31.2 



Ethanol 

.50 

1.0 



46. 1 

1516 

37.4 



Freon TF 

.50 

1.0 



31.9 

3490 

43.2 



Ethanol 

.75 

1 

.5 



41.2 

2036 

44.7 

I 

Ethanol 

.50 

1.0 



52.3 

1720 

48.2 

T 

Freon TF 

.50 

1.0 



37. 5 

4230 

59.5 

I 


Ethanol 

.75 

1 

5 



47.5 

2347 

59.6 

I 




.50 

1.0 

1 


85.7 

2820 

129.0 

I 




.75 



.75 

12.7 

630 

4.2 

S 









22.5 

1110 

13.3 

S 









24.8 

1227 

16.2 

£ 










27.3 

1350 

19.7 


r 









29. 4 

1454 

22.9 











32. 0 

1580 

27.0 











35.0 

1730 

32.3 











39.1 

1930 

40.2 

1 



a Do not satisfy minimum -Reynolds- number criterion, but included as data 
and in data points. 




Figure 1. - Schematic of liquid jet impinging on a flat plate. 



Figure 2. - Flow pattern as a function of nozzle 
height. 



Figure 3. - Flow pattern as a function of jet radius. 
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Figure 8. - Schematic diagram 
of sharp-edge disk. 



(a) Surface tension flow. 


(b) Transition flow. 


(c) Inertial flow. 


Figure 9. - Liquid jet impinging on flat solid disk. 
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Figure 10. - Zero-gravity experimental results. 
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(a I Liquid jet impingement. 
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(b) Surface tension. 

Figure 11. - Models in physical plane. 
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Figure 12. - Circular liquid jet impinging 
on infinite plate. 
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(b)z Solution. 

Figure 13. - Inverse formulation excluding surface tension for finite plate. 
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Figure 15. - Nodal point representation for rec 
tangular mesh. 
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(b)z Solution. 

Figure 14. - Surface tension model - inverse formulation for finite plate. 












(a) r Solution. 

Figure 17. - Finite difference representation for surface-tension -dominated model, 













Dimensionless axial coordinate, z/Rq 



Dimensionless radial coordinate, r / R 0 ment on finite plate, excluding surface tension, for ratio of 

jet radius to disk radius Rn/L of 1/2. 

Figure 20. - Numerical solution of liquid jet impinging on finite plate, ” 

excluding surface tension, for ratio of jet radius to disk radius Rq/L 
of 1/2. 






Figure 24. - Comparison of numerical results for infinite plate with ref- 
erence 29. 



Figure 25. - Comparison of numerical results for finite- 
plate inertial flow with experiments for ratio of jet ra- 
dius to disk radius Rq/L of II Z 
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Figure 26. - Comparison of numerical results 
for finite-plate inertial flow with experiments 
for ratio of jet radius to disk radius Rg/L of 
3/4. 
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0 .25 .50 .75 

Ratio of jet radius to disk radius, Rg/ L 

Figure 27. - Experimentally obtained free -surface data 
showing downward trend of jet at higher ratios of jet 
radius to disk radius. Dimensionless radial coordin- 
ate, r / Rg, 4. 





/“ Wire-release mechanism 
Anvil -v / double-acting piston 
\ / access hole 



(a) Before release. (b) During free fall. (c) After deceleration. 

CD-8418 

Figure 29. - Position of experiment package and drag shield before, during, and after test drop. 
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